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The  purpose  of  study  was  to  investigate  the  application 
of  theories  for  the  design  of  robust  control  systems  for 
multiple  input  multiple  output  systems,  and  to  apply  these 
techniques  to  a  realistic  control  system.  The  techniques 
investigated  are  known  as  the  and  design  techniques,  and 
involve  the  minimization  of  the  2-  and  oo-norms  of  the  system 
transfer  function  matrices.  Considerable  system  modification 
is  required  to  apply  these  techniques,  and  a  methodology  for 
accomplishing  this  modification  is  developed.  Evaluation  of 
the  controllers  developed  requires  simulating  the  system 
performance,  and  a  digital  simulation  is  developed  for  this 
purpose . 
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and  design  techniques  are  applied  to  a  state  space 
representation  of  a  surface-to-air  missile  system  with 
unstable  nonminimum  phase  airframe  response  characteristics. 

A  method  for  converting  the  tracking  and  command  following 
problems  to  the  "standard  problem"  is  given.  Artificial 
frequency-dependent  filtering  of  external  inputs  is  required 
for  proper  closed  loop  system  performance,  and  a  method  for 
deriving  the  filters  is  developed.  Areas  requiring  further 
study  are  identified. 


DESIGN  OF  A  ROBUST  CONTROLLER  FOR  AN  UNSTABLE 
NONMINIMUM  PHASE  GUIDED  MISSILE 

L-  INTRODUCTION 

The  purpose  of  this  thesis  is  to  investigate  and  apply 
new  and  robust  control  system  design  techniques  to  a 
surface-to-air  missile  control  system,  and  compare  the 
performance  of  the  resultant  systems  to  the  performance  of  a 
system  with  a  controller  derived  from  conventional  control 
system  synthesis  techniques. 

A  surface-to-air  missile  example  was  chosen  for 
comparison  because  it  will  test  the  robust  controller  design 
methodologies'  capability  tc  handle  extreme  cases.  The 
transfer  functions  of  the  chosen  missile  have  right  half-plane 
poles  and  zeros  (open-loop  unstable  and  nonminimum  phase). 

The  missile  is  unstable  due  to  the  location  of  the  center  of 
gravity,  which  is  aft  of  the  aerodynamic  center,  and  is 
nonminimum  phase  due  to  the  use  of  tail-mounted  control 
surfaces.  The  objective  is  to  design  a  controller  which 
stabilizes  the  system  and  makes  it  responsive  to  commands  so 
that  guidance  commands  successfully  steer  the  missile  to  the 
target . 

Three  design  methodologies  will  be  compared.  The  first 
is  a  classical  approach,  using  root  locus  and  time  response 
analyses  to  derive  a  control  system  with  good  performance  in 
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the  .-.•-■minai  case.  The  second  is  the  H,  approach  which 

provides  robustness  for  systems  exposed  to  well-defined 

stochastic  disturbances  and  parameter  variations.  The  third 

is  t-.He  H  methodology  which  minimizes  system  response  to 
00 

unknown  but  deterministic  disturbances  and  provides  robustness 
with  respect  to  unstructured  uncertainties.  The  relative 
performance  of  these  methods  will  be  compared  by  implementing 
each  controller  in  a  two  degr ee -o f - freedom  digital  simulation. 
To  provide  a  baseline  for  performance,  each  controller  will 
first  be  implemented  in  a  nominal  model  which  is  not  subject 
to  disturbances.  Each  will  then  be  implemented  in  a  more 
realistic  model  which  provides:  (a)  both  time-varying  and 
constant  parameter  variations  in  the  system  components  and 
equations  of  motion,  (b)  stochastic  disturbances  via  target 
glint,  sensor  noise  and  turbulence,  and  (c)  deterministic 
disturbances  via  loss  of  target  data.  The  measure  of  merit 
used  for  these  simulation  runs  will  be  miss  distance,  as  well 
as  detailed  examination  of  missile  transient  responses. 
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ii_.  Overview 

The  scope  of  this  report  is  broad,  as  a  large  number  of 
tasks  must  be  completed  prior  to  comparing  the  performance  of 
the  different  controllers.  These  include  development  of  the 
system  description,  conversion  to  state  space  format, 
derivation  of  the  different  controllers,  and  simulation  of  the 
missile  engagement  scenarios.  This  chapter  details  how  these 
tasks  are  organized  and  presented  in  this  report,  so  it  will 
be  easier  for  the  reader  to  follow  the  development. 

Mi sal Is  System  Description 

In  Chapter  III,  the  missile  system  is  described  in 
detail.  The  equations  of  motion  for  the  missile  body  are 
developed  and  the  system  transfer  functions  are  given  for  both 
the  missile  body  and  the  target  tracking  system.  Flight 
conditions  are  described  for  both:  (1)  a  nominal  situation 
where  the  environment  and  missile  system  are  not  subject  to 
changes  or  variations,  and  (2)  a  more  realistic  situation 
where  the  environment  changes,  the  missile  system  components 
are  subject  to  variation,  and  disturbances  enter  the  system. 

Space  Representation 

In  Chapter  IV,  the  missile  system  description  given  In 
Chapter  III  is  converted  to  a  state  space  representation, 
which  is  necessary  to  apply  the  robust  control  methods.  Both 
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the  nominal  model  and  the  model  with  variations  are  described 

in  detail. 

Classical  Control  System  Design 

In  Chapter  V,  a  controller  for  the  missile  system  is 
derived  using  standard  classical  design  techniques,  including 
root  locus,  time  response,  and  frequency  domain  analysis. 
Acceptable  performance  is  obtained  with  a  simple  proportional 
feedback  controller.  The  author  realizes  that  better 
performance  could  probably  be  obtained  with  lead/lag, 
integrator,  and  derivative  controllers.  Due  to  scope 
limitations,  however,  these  options  were  not  pursued  in  favor 
of  focusing  on  robust  controller  design. 

H2  Controller  Design 

In  Chapter  VI,  the  theory  for  H2  controller  design  is 
reviewed  and  applied  to  the  missile  system.  This  involves 
some  manipulation  of  the  original  state  space  representation 
and  addition  of  many  artificial  weighting  filters  to  the  plant 
Inputs  and  outputs  to  get  the  problem  into  the  form  solvable 
by  existing  theory.  Producing  these  manipulations  and 
artificial  filters  constituted  the  bulk  of  the  design  effort, 
and  the  methods  for  selecting  these  are  described  in  detail. 

a,*  contrail  sr  Efi&lsa 

In  Chapter  VII,  the  theory  for  controller  design  is 
applied  to  the  missile  system.  The  plant  manipulations  and 
artificial  filtering  of  Inputs  and  outputs  needed  for  this 
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design  effort  are  similar  to  those  needed  for  the  design. 
For  this  system,  no  additional  problem  formulation  effort  is 
required,  so  this  chapter  is  devoted  solely  to  a  review  of  the 
problem  requirements  and  derivation  of  the  controller. 

Performance  CompaxlS-Qiia. 

In  Chapter  VIII,  the  performance  of  each  of  the 
controllers  developed  is  compared.  Comparison  is  based  on 
missile  flyout  simulation  results,  with  miss  distance  being 
the  primary  figure  of  merit.  Development  of  the  flyout 
simulation  is  briefly  discussed.  Comparison  of  flyout  time 
histories  allows  some  qualitative  comparisons  in  addition  to 
the  quantitative  results  yielded  by  miss  distance  comparison. 

Results  and  Conclusions 

In  Chapter  IX,  the  author  makes  remarks  about  the 
advantages  and  disadvantages  of  the  robust  controller  designs 
and  methodologies.  Areas  of  further  study  are  identified. 

Filters 

The  artificial  filters  developed  for  the  and 
controller  designs  are  described  in  detail  in  Appendix  A. 

Flyout  Time  Hlatorlsa 

Appendix  B  contains  missile  system  output  time  histories, 
missile  and  target  position  plots,  and  tabulated  state 
histories  for  a  representative  sample  of  the  flyout  runs. 
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The  controllers  obtained  for  the  H_  and  H  designs  are 
composed  of  24  states  each.  Reduced-order  controllers  were 
obtained  for  use  in  the  simulation.  Both  the  full-order  and 
reduced-order  controllers  are  described  in  detail  in 
Appendix  C. 


Ill-  £  UI  t  '1C  c  ~  1 0  -  A 1 1  Mluwl  lc  System  Description 

The  surface-to-air  missile  system  description  used  herein 
is  provided  by  Sverdrup  Technology,  Inc.  [12].  The  system 
description  is  derived  from  a  digital  model  generated  by  Jim 
Coburn  of  Sverdrup  as  a  generic  model  of  a  tail-control 
missile  to  be  used  for  training  and  experimentation  purposes. 

The  model  is  composed  of:  (1)  a  system-specific  data 
base,  and  (2)  a  simulation  execution  program.  The  simulation 
execution  program  is  generic  in  its  routines,  and  depends  upon 
specific  data  provided  by  the  data  base  for  simulation  of  a 
particular  missile  system.  The  data  base  for  a  generic 
tail-control  missile  created  by  Coburn,  with  some  slight 
modifications,  was  used.  These  modifications  include 
simplification  of  the  antenna  servomechanism  for  ease  of 
computation,  omission  of  lovpass  filters  associated  with  the 
servo  op-amp,  moving  the  center  of  gravity  aft  to  cause  the 
system  to  be  statically  unstable,  and  simplification  to  two 
dimensions . 

Missile  velocity  and  altitude  time  histories  are 
postulated  to  provide  3ome  realistic  parameter  variations. 

Any  Inaccuracy  in  this  representation  of  the  model  is  solely 
the  fault  of  the  author. 

The  model  developed  herein  describes  motion  in  a  single 
plane  only,  and  does  not  simulate  missile  axial  forces  (drag) 
or  the  effects  of  gravity. 


3-1 


Equation  a  oJL  Motion 

The  following  equations  of  motion  [2]  model  the  airframe 
dynamics.  They  are  used  to  derive  the  transfer  functions  and 
state  space  models  associated  with  the  airframe. 


V3’ 

=  N  a(  s  ) 
o 

+  N6  6( 5) 

(3.1) 

V31 

=  -Ma  a(s) 

+  <5(s) 
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<*(s ) 

=  Vs)  - 

y(s) 

(3.3) 

r(  s ) 

.  vs)  . 

N  a(s  )  -t-N ,  <5(s  ) 

a  o 
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m  iv| 

with  the  above  defined  as 

=  component  of  acceleration  normal  to  airframe 
a  missile  heading  angle 
«  *  angle  of  attack 

r  -  bearing  of  missile  velocity  vector 
6  =  control  surface  deflection 

| V |  =  magnitude  of  missile  velocity 

The  dimensional  stability  derivatives  [9:412-4451  are  given  by 


N  =  C 
a  n 

q 

*re£/m 

(3.5) 

*  C»I 

q 

‘ref7" 

(3.6) 

M  =  C 

a  m 

a 

q 

»ref  Lref71 

(3.7) 

M  &  =  C 

6  n<5 

q 

*re£  Lref 71 

(3.8) 

where  the  stability  derivatives  (Cn  ,  etc.)  are  tabulated  in 

a 

[121  as  a  function  of  Mach  number,  and  the  other  parameters 
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it  ^  C*  6  Z  i  n e  d  U 3 

2  2 

q  =  dynamic  pressure  ( =pV  /2)  (N/m  ) 

2 

A  ,  =  missile  reference  area  (=  0.031m  t 1 2 ] ) 
ref 

m  =  missile  mass  (kg) 

2 

I  =  missile  moment  of  Inertia  (kg-m  ) 

Lref  =  Missile  reference  length  (=  0.2m  [12]) 

The  sign  conventions  used  for  the  missile  airframe  parameters 
are  as  shown  in  Figure  3.1. 


The  target  tracking  antenna  and  servomechanism  are 
attached  to  the  front  end  of  the  missile.  The  tracker  data  is 
given  in  transfer  function  form  in  C 1 2 1  and  will  be  listed  in 
the  next  section. 

Nominal  Flight  condition 

This  section  lists  the  values  used  in  the  nominal  model. 
For  this  time-invariant  description,  all  flight  conditions  are 
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held  cunstant,  and  no  noise  or  disturbances  enter  the  system. 

The  nominal  flight  conditions  (fixed  based  on  data  from  (12] 
and  with  atmospheric  data  from  (11)  are  listed  in  Table  3.1. 

A  computer  simulation  of  the  missile  flyout  will  be 
constructed  using  these  constant  values  for  the  nominal  flignt 
condition.  This  model  will  be  used  as  a  baseline  for 
determining  noise  and  disturbance  effects. 

Table  3.1.  Nominal  Flight  Conditions 


Center  of  Gravity 

2.63  m 

Mass 

127  kg 

Inertia 

271.3  kg-.m 

Altitude 

20  kft 

Mach 

3.5 

Velocity 

1100  m/s 

Speed  of  sound 

315  m/s 

Dynamic  pressure 

4xl05  N/m2 

Air  density 

0.65  kg/m3 

Nominal  Airframe .  The  airframe  stability  derivatives  for 
thi3  nominal  flight  condition  were  extracted  from  the  tabulated  data 
in  [12]  and  are  given  in  Table  3.2. 


Table 


Stability  Derivatives 

11.5 

5.7 

-0.1 

-0.3 


The  values  in  Tables  3.1  and  3.2  were  then  applied  to 
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equations  3. 5-3. 3  to  yield  the  dimensional  stability 
derivatives  Hated  in  Table  3.3. 


Table  3.3.  Dimensional  Stability  Derivatives 

2.0  (g's/deg) 

N.  1.0  (g's/deg) 

°  2 

-0.91  (deg/sec  /deg) 

-2.74  (deg/sec^/deg ) 


These  coefficients  yield  the  following  transfer  functions, 
using  the  equations  of  motion  listed  above  and  adding  the 
servomechanism  dynamics 


<5(s) 


3  (s) 

com 

(S) 

m 

6  (s ) 
com 

0  (S) 
m 

3  (s ) 
com 

q(s ) 

3  (s) 

com 


Vs) 

6com(sr 


77 

s  +  77 


(12] 


-211(3+1.1894) 

"  (s+77) (s+1. 5917) (s-0. 5717) 

-211(3+1.1894) 
s( s+77) (s+1.5917) (s-0.5717) 

-39.3(s+5.3725) 

(s+77) (s+1.5917) (s-0.5717) 

77(3+2.5278) (s-2.5278) 
(3+77) (s+1.5917) (s-0.5717) 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


HLaCklna  S^aisa.  The  target  tracking  system  will 


be  synthesized  from  two  components;  a  simplified  boresight 
error  resolver  and  an  antenna  mounted  on  a  servomechanism. 

The  boresight  error  resolver  returns  an  error  signal,  4>,  in 
response  to  the  boresight  error  input  <pQ  (angular  position  of 
target  relative  to  antenna  boresight  position  ),  with  a 
transfer  function 
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P  (3)  =  20.7 

0e(s)  “  s+5.26 

The  sign  convention  for  <p  is  as  shown  in  Figure  3.2 
is  the  angular  displacement  of  the  antenna  relative 
missile  centerline  and  <p^  is  the  bearing  to  target. 


(3.14) 

,  where  & 

'  s 

to  the 


The  antenna  servomechanism  Is  attached  to  the  antenna 
pedestal  such  that  It  generates  an  angular  displacement  rate 
of  the  antenna  boresight  0  in  response  to  an  input  signal  «. 

9 

The  transfer  function  for  the  servomechanism  is 


V8>  =  94 

w(s)  s+104 


(3.15) 


This  combined  servoaotor/resolver  system  can  be  represented  as 
shown  in  Figure  3.3. 
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*Tls) 


imvmm 

Vs) . 

i 

Vs>  Ti+  Vs) 

20 . 7 

0(S) 

IKrnia 

s 

f°  I 

s+5. 26 

Figure  3.3.  Target  Tracking  System 


Coupled  System.  When  the  antenna  is  coupled  to  the 
airframe  the  total  integrated  system  looks  as  in  Figure  3.4 


<*>T(s ) 


com 


is) 


m 


Figure  3.4.  Coupled  System  Block  Diagram 


Note  that  the  output  of  the  error  resolver  has  not  yet  been 
fed  into  the  antenna  servomechanism.  This  will  be  done  during 
the  design  process. 

Derivation  of  the  additional  transfer  functions 
associated  with  the  coupled  system  yields 


3- 


P\  3  '/ 


4353(5+1.1334) 


(  3  .  x  6  j 


6  (s)  s(s+77) (s  +  1. 5917) (s-0. 5717) (s+5. 26) 

com 

4>{s  )  =  20 . 7 

4>t(s)  (s  +  5. 26) 

<p(s  )  _  -1945.8 _ 

w(s)  s ( s+104 ) ( s+5 . 26 ) 

Time.  Varying  System  yit h  Disturbance  Inputs 

This  section  lists  all  the  time-varying  parameters, 
parameter  variations,  and  disturbance  inputs  to  the  system. 
These  effects  are  simulated  in  the  time-varying  model. 

Missile .  In  reality,  the  missile  does  not  operate  in  a 
constant  deterministic  environment  and  many  of  the  physical 
parameters  of  the  system  are  not  known  accurately  due  to 
component  variation  from  missile  to  missile.  A  realistic 
simulation  of  the  missile  mi’3t  include  the  effects  of  these 
variations.  To  this  end,  time  functions  will  be  used  to  model 
parameters  vhich  vary  over  the  time  of  guided  operation  (8 
seconds),  some  of  the  component  properties  will  be  perturbed, 
and  disturbances  will  be  added. 

Time  Varying  Parameters .  In  order  to  realistically 
vary  flight  conditions  with  time,  the  assumption  is  made  that 
the  missile  has  finished  boosting,  initiated  the  sustained 
propulsion  phase,  and  starts  guidance  at  a  15kft  altitude, 
climbing  at  a  constant  rate  of  1000  feet  per  second.  This 
flight  path  is  assumed  based  on  a  typical  surface-to-air 
missile  launch,  and  the  target  altitude  induced  by  this 


(3.17) 

(3.18) 
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profile  is  dependant  upon  actual  time  to  intercept.  The  time 
function  for  the  missile  altitude  is  then 

h ( t )  =  15,000  +  lOOOt  (feet)  (3.19) 

Given  standard  atmospheric  conditions,  this  causes  the  air 
density  and  local  speed  of  sound  to  vary  as  [1:701 

p(t)  =  0.7626  -  0 . 0214t  (kg/m3)  (3.20) 

a(t)  =  322.27  -  1.26t  (m/s)  (3.21) 

Missile  velocity  normally  varies  slightly  as  the  propellant  is 
burned  during  the  sustained  phase,  so  the  following  velocity 
time  function  is  assumed  based  on  propulsion  data  from  [121 

V(t)  =  1000  +  80t  -  8t 2  (m/s)  (3.22) 


Since  the  reference  data  tabulates  aerodynamic  coefficients  as 
a  function  of  local  Mach  number,  it  is  necessary  to  express 
this  value  as 


M(t)  =  V(t)/a(t )  = 


1000 

322.: 


80t  -  8t' 
r  -  1.26t 


(3.23) 


For  ease  of  calculation,  it  is  desired  to  express  this 
function  as  a  simple  polynomial  in  time  only.  This  can  be 
accomplished  by  using  a  Taylor  series  expansion  (which  for  the 
case  of  ratios  of  polynomials  is  equivalent  to  carrying  out 
the  division).  The  Taylor  series  expansion  yields  an  infinite 
order  polynomial,  which  is  normally  truncated  at  the  point 
where  the  desired  accuracy  is  reached  for  the  interval  of 
interest.  An  interval  of  eight  seconds  and  an  accuracy  of  2% 
are  used  for  all  of  the  Taylor  series  expansions  developed  in 
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this  section  (that  is,  when  the  highest  order  term  contributes 

less  than  2%  to  the  value  at  eight  seconds,  it  and  all 
following  terms  are  ignored).  Taylor  series  expansion  yields 
the  Mach  number  time  function 

M ( t )  =  3.103  +  Q.26t  -  0.0238t2  (dimensionless)  (3.24) 


Based  on  the  above  values  for  density  and  velocity,  dynamic 
pressure  can  be  expressed  as 


q ( t )  =  p( t)  V ( t ) 2/2 

=  (0.7626  -  0.0214t) <1000  +  80t  -  8t2)2/2  (3.25) 


Inertia,  mass,  and  stability  derivatives  were  derived  from  the 
tabulated  data  [12],  where  they  are  listed  as  function  of  Mach 
number  and  time.  Using  equation  3.24,  these  were  all 
converted  to  time  functions,  yielding 


m(t) 

=  127  -  2t  (kg) 

Acg ( t ) 

=  0.84  -  0 . 0008t 

(m) 

I  ( t ) 

=  279.6  -  1 . 58t  + 

0 . 0027t2 ( kg-m2 ) 

C  (t) 

=  11.30  -  0 . 0447t 

+  0.0409t2 

n 

a 

C  (t) 

=  6.130  -  0 . 0447t 

+  0 . 0409t2 

C  (t)  =  0.0017  -  0 . 037t  +  0 . 0038t2 

at 

a 

C  (t)  =  -0.329  +  0 . 0181t  -  0 . 0017t2 
m<5 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


Using  equations  3. 5-3. 8  for  dimensional  stability  derivatives, 
equations  3.25-3.32  are  used  to  generate  the  following 
equations  for  dlmensionsal  stability  derivatives  as  functions 
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of  time 


N  (t)  =  1.871  +  0 . 202t  -  0.0262t2+  0.0008t3  (3.33) 

a 

N  ,  ( t )  =  1.016  +  0 . 076t  -  0 . 0161t2  -  0.0004t3  (3.34) 

o 

M  (t)  =  +0.0144  -  0 . 3109t  -  O.Olllt2  +  0.0085t3  (3.35) 

ot 

M6(t)  =  -2.78  -  0 . 230t  +  0.044t2  -  O.OOlt3  (3.36) 

Component  Var latlons ■  Many  other  parameters  of  the 
system  vary  from  missile  to  missile,  but  not  with  time.  In 
composing  these  variations,  the  author's  best  estimate  was 
used . 

The  time  constant  (rf)  of  the  fin  servomechanism  can 
vary.  The  perturbed  servomechanism  transfer  function  can  be 
written  as 

<5(3)  =  1/Tf  =  77  +A  (3.37) 

<5  m(s)  s  +  1/T  s  +  ( 77+A) 
com  r 

where  A  is  normally  distributed  with  zero  mean  and  a  variance 
2 

a  determined  by  manufacturing  quality  control.  A  variance 
value  of  4  was  assumed.  Standard  statistical  notation  denotes 
this  as  A  =  N ( 0 , 4 ) .  Similarly,  it  was  assumed  that  the  body 
rate  sensor  suffers  from  a  variance  of  1%  in  gain,  so  it  sees 
a  variation  of  N(0,.01)  from  Its  normal  gain  of  1.  Also,  It 
was  assumed  the  normal  acceleration  sensor  sees  a  variation  of 
N ( 0, . 01 )  from  its  normal  gain  of  1. 

Disturbances .  There  are  disturbances  in  the  form  of 
turbulence  and  power  supply  noise  entering  the  missile 
airframe  system.  Reference  19)  presents  data  describing 
turbulence  disturbances,  and  from  this  reference  atmospheric 
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turbulence  can  be  modeled  aa  random  wind  -gusts  of  velocity  w 

d 

where  wg  =  N(0/3)  (m/sec).  Also  from  [9],  it  can  be  seen  that 
wind  gusts  can  be  modeled  as  affecting  the  airframe  as  a 
change  in  angle  of  incidence.  Specifically,  the  wind  gust 
can  be  modeled  as  an  additional  component  to  a,  namely  a  , 
with  the  transfer  function 

a  ( s  )  . 

v*(s)  =  | V | ( s  +  1/t  )  (3.38) 

g  11  m 

where  is  a  time  constant  dependant  on  the  airframe.  It  is 
assumed  for  this  missile  that  is  very  large,  such  that  the 
above  transfer  function  can  be  approximated  a3 


Vs)  .  1 

V  (s)  '  |V|s 


ag(s) 


w  (s) 

1 W~ 


(3.39) 


(3.40) 


Using  the  nominal  velocity  as  an  approximation,  the  turbulence 
disturbance  is  modeled  as  a  =  N(0, 0.003). 

g 

Another  disturbance  entering  the  missile  system  is  power 

supply  noise  (?).  This  disturbance  appears  as  a  random  input 

to  the  fin  servomotor,  entering  at  the  same  point  as  <5com. 

Such  noises  can  be  characterized  as  having  zero  mean  with  some 
2 

variance  a  dependent  on  the  magnitude  of  the  noise  in  the 

system.  This  noise  also  has  distinct  frequency 

characteristics,  and  can  be  represented  as  a  random  signal  ft  - 
2 

N ( 0 ,9  )  being  passed  thru  a  filter  with  the  same  frequency 


3-12 


character  is t ics  as  the  noise.  For  instance,  if  the  noise  lies 

within  the  frequency  range  1 ,  it  can  be  represented  as 


?(s) 


s  ( s  +  )  ( s  +  o>2  ) 


(3.41) 


where  Kf  is  set  to  achieve  unity  gain  over  the  frequency  range 
of  interest.  For  the  power  supply  noise,  the  following  values 
are  assumed 


ft  =  N ( 0, 0 . 05 )  (3.42) 

oo^  =  oo^  =  400  sec  ^  (3.43) 


Target  Tracking  SysJtfiJtt.  The  target  tracking  system, 
being  composed  of  electronic  parts  and  not  particularly 
dependent  upon  the  flight  conditions  for  its  operating 
parameters,  does  not  suffer  any  of  the  time-varying  parameter 
changes  that  the  missile  does.  It  does  see  component 
variations  and  disturbances  similar  to  the  missile,  however, 
and  these  must  be  addressed. 


The  resolver  can  be  modeled 


as  consisting  of  a  time  constant  r r  and  a  gain  K^, 

#(s)  _  r 

*e(S)  "  S  +  1/T^ 


(3.44) 


and  is  subject  to  variations  in  both  of  these  values.  Kr  is 
assumed  to  vary  as  N(0,2)  and  1/Tr  as  N(0,0.5).  The 
servomotor  can  be  modeled  in  this  fashion  also.  For  the 


servomotor,  it  is  assumed  that  Kg  varies  as  N(0,4)  and  1/Tg  as 
N ( 0 , 4 )  (Independently).  In  addition,  it  was  assumed  that  the 
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censors  which  detect  the  rate  and  position  of  the  antenna 
suffer  variances  of  1%  in  accuracy,  so  each  of  these  sees  a 
variation  of  N(0,0.01)  from  its  normal  gain  of  1. 

Disturbances .  The  electrical  components  of  the 
target  tracking  system  are  subjected  to  the  same  power  supply 
noise  as  the  fin  servomotor  described  above.  Thus,  this  same 
noise  will  be  injected  into  the  angle  error  resolver  and  the 
antenna  servomotor. 

Finally,  target  data  is  corrupted  by  target  glint.  Glint 
enters  the  system  as  a  random  pulse-to-pulse  variation  in 
target  position.  The  glint  was  assumed  to  be  distributed  as 
N ( 0 , 2 5 )  (meters).  This  can  be  converted  to  an  angular  noise 
input  £  dependent  on  target  range  (R^,)  by  using  the  small 
angle  approximation  and  converting  radians  to  degrees  to  get 

K  =  N(0,1500/Rt)  (3.45) 

All  of  the  above  noise  and  disturbance  inputs  are  modeled  in 
the  simulation.  To  do  this,  it  Is  necessary  to  model  how 
these  inputs  enter  the  system  in  state  space  form.  This  is 
done  in  Chapter  IV. 
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The  state  space  representation  of  the  airframe  is  derived 
from  the  equations  of  motion.  T^e  standard  state  space 
nomenclature  [6:21 

x  =  A  x  +  B  u  (4.1) 

y  =  C  x  +  D  u  (4.2) 

will  be  used  as  a  base  for  building  the  state  space 
representations  used  herein.  The  systems  described  will  grow 
too  large  to  maintain  this  simple  notation  and  preserve 
clarity,  so  subscripts  and  additional  vectors  will  be  added  as 
needed.  The  development  of  this  representation  will  parallel 
the  Laplace  domain  development  in  Chapter  III  as  much  as 
possible . 


Missile.  Aisisaae. 

Using  the  standard  nomenclature  above,  the  missile 
airframe  is  represented  by 


x  = 


<5 

e 

■ 

e 

a 

a 


(fin  deflection,  deg) 
(body  rate,  deg/sec) 
(heading  angle,  deg 
(angle  of  attack,  deg) 


(4.3) 


u  *  6  (fin  deflection  command,  deg)  (4.4) 

coa 


y  3 


(normal  acceleration,  g's) 


(4.5) 
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Using  these  state,  input,  and  output  vectors,  the 
corresponding  A,  B,  C,  and  o  matrices  from  the  transfer 
functions  in  Chapter  III  are 


-77  0  0  0 

-2.74  0  0  0.91 

0  10  0 
-0.51  1  0  -1.02 


(4.6) 


B  = 


(4.7) 


10  0  2 
0  10  0 


(4.8) 


D  = 


0 

0 


(4.9) 


This  system  is  controllable  [6:531,  but  the  state  is 
not  observable  [6:661. 

Tat-ast.  Tracking  System 

For  the  target  tracking  system,  the  vectors  are 


u  = 


Oi 


(boresight  error  signal,  deg) 

(antenna  rate,  deg/sec) 

(antenna  angle  measured  from 
missile  centerline,  deg) 


(4.10) 


(bearing  to  target,  deg) 
(antenna  servo  input,  deg/sec) 
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(4.12) 


(4.13) 


(4.14) 


(4x5) 


(4.16) 


Coupled  SysLtea 

Combining  the  tracker  with  the  airframe,  and  making  the 
substitution  ^e=^~0m~03  (See  Figure  3.4)  yields  the  following 
state-space  representation  of  the  combined  system,  hereafter 
referred  to  as  the  plant.  At  this  time,  it  is  convenient  to 
split  the  input  vector  into  two  separate  vectors,  one  of  which 
contains  those  inputs  which  can  be  manipulated  by  the 


designer,  &com  an<*  the  other  containing  the  input  which  is 
an  external  signal,  These  input  vectors  will  hereafter  be 

designated  u  and  w,  respectively.  Note  the  additional 
subscripts:  p  to  denote  the  plant,  11  to  denote  those  matrices 


corresponding  to  the  input  u 


P' 


and  12  to  denote  those  matrices 
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corresponding  Co  Che  input  w  .  Thus,  the  plant  i 

by 

x  =  A  X  +  B  ..  ..  u  +  B._w 
P  P  P  PH  P  pl?  P 

y  =  C  x  +  D  ,  u  +  D._v 
JP  P  p  pH  P  pl2  p 

where 


x 


P 


u 


-77 

-2.74 

0 

-0.51 

0 

0 

0 


6 

e 

m 

9 

m 

a 

<p 

e 

s 

9 

s 


<5 

com 

U) 


[  <Prr  \ 


m 

An 

<P 

9 

s 


0  0  0.91  0  0  0 

1  0  0  0  0  0 

1  0  -1.02  000 

0  -20.7  0  -5.26  0  -20.7 
0  0  0  0  -104  0 

0  0  0  0  1  0 


s  represented 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 
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(4.24) 


(4.25) 


(4.26) 


(4.27) 


0 
0 

0  (4.28) 

0 
0 

% 

This  is  the  'clean'  version  of  the  missile  model  and  will 
be  used  for  the  classical  design  approach. 


Time.  VaxYlJDfl  Syatea  with 


Inputs 


In  addition  to  the  above  system  representation,  the 
realistic  system  is  subject  to  time-variant  changes  in 


environmental  conditions,  component  variations,  and  noise 
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inputs  as  described  in  Chapter  III. 

Disturbances .  The  noise  inputs  are  appended  to  the  clean 
state  space  representation  such  that 


=  A  X 

P  P  P 

*  bpu 

u 

P 

+  Bpl2  "p  + 

B  - 
p2 

* 

(4.29) 

v  =  C  x 

+  D  .  , 

u 

+  D  .  «  w  + 

D  . 

£  +  T) 

(4.30) 

P  P 

pll 

p 

pl2  p 

p2 

where  equations 

4.17-4 

.22 

still  apply 

and 

B 


P2 


P2 


0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 


(4.31) 


(4.32) 


?  = 


V  - 


N (0,0. 003) 

N ( 0, 1500/R- ) 
N ( 0, 0 . 05 )  1 


N  (  0 , 0 . 1 ) 

N ( 0, 0 . 05 ) 
0 

N(0, .1) 

N(0, .05) 


( turbulence ) 

(target  glint) 

(power  supply  noi3e) 

(body  rate  gyro  noise) 
(accelerometer  noise) 
(no  measurement  noise) 
(rate  gyro  noise) 
(potentiometer  noise) 


(4.33) 


(4.34) 


Two  items  in  the  above  deserve  note.  First,  notice  that 
target  glint  has  a  variance  that  is  a  function  of  range. 

This  Indicates  that  this  signal  is  non-stationary,  which 
creates  special  problems.  This  is  another  time-varying 
parameter.  Secondly,  notice  that  t?'s  third  element  is 
zero.  This  is  because  <P  is  an  electrical  signal  directly 
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available  for  measurement,  and  has  no  measurement  no^se 
corruption  associated  with  it. 

In  addition  to  these  random  noise  inputs,  the  real  system 
is  also  subjected  to  a  discrete  deterministic  (but  unknown) 
disturbance,  specifically  loss  of  target  data  due  to  ECM. 

While  qualitative  statements  can  usually  be  made  about  factors 
such  as  this,  it  cannot  be  treated  as  a  known  quantity  which 
can  be  compensated  for. 

Time-varying  parameters  and  Component  Variations .  Next, 
parameter  variations  must  be  considered.  These  will  enter 
into  the  real  system  as  variations  in  missile  velocity,  mass, 
inertia,  atmospheric  conditions,  center  of  gravity  location, 
servomechanism  time  constants,  component  gains,  and 
high-frequency  unmodeled  dynamics.  These  effects  will  all  be 
simulated,  except  the  high-frequency  unmodeled  dynamics.  This 
exception  is  due  solely  to  the  desire  to  limit  the  scope  of 
this  investigation.  This  variation  of  parameters,  both 
time-varying  and  invariant,  will  be  modeled  by  making  the  A, 

B,  and  C  matrices  functions  of  time  and  adding  the  component 
variations  as  follows 


-77+ A 

0 

0 

0 

0 

0 

0 

0 

0 

-M 

a 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-H6/v 

1 

0 

-Na/V 

0 

0 

0 

(4.35) 

0 

0 

-20 .7+A2 

0  -5. 

26  +  A3 

0  - 

20.7+A2 

0 

0 

0 

0 

0 

104+A5 

0 

0 

0 

0 

0 

0 

1 

0 
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7  7*A  0 

0  0 

0  0 

0  0 

0  0 

0  9  4  +A 

4 

0  0 

L_  J 


(4.36) 


B 


pl2 


0 

0 

0 

0 

20.7  +-^2 
0 
0 


(4.37) 


0  i+a7 
N<5(1+V  0 

Cp  -  0  0 

0  0 

0  0 


0  0  0 

0  N  (1+AC)  0 

oi  b 

0  0  1 

0  0  0 

0  0  0 


0  0 

0  0 

0  0 

1+A0  0 

08  1*A, 


(4.38) 


where  N^,  N^,  M^,  and  V  vary  as  defined  in  Chapter  III  and 

the  A's  are  defined  as 


A 


N ( 0, 4  ) 
N(0,2) 

N ( 0, 0 . 5 ) 
N<0,4) 

N ( 0, 4  ) 
N(0,0.01) 
N( 0, 0 . 01 ) 
N( 0, 0 . 01 ) 
N ( 0, 0 . 01 ) 


(4.39) 


Note  that  the  A's  are  preset  and  do  not  vary  with  time. 
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With  the  problem  now  in  the  proper  form,  the  design 


effort  Is  initiated  in  Chapter  V. 
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'£•  Zi&sS-i&AL  v  v ill  £  u  1  jys  Design 

In  this  section,  the  objective  is  to  derive  a  missile 
control  system  based  on  classical  design  principles.  Using 
the  root  locus,  bode  plot,  and  step  response  analysis  tools,  a 
controller  is  designed  which  stabilizes  the  missile,  provides 
good  tracking  loop  performance  for  the  target  tracker,  and 
guides  the  missile  to  the  target  location  using  standard 
proportional  navigation  techniques. 

The  standard  methodology  for  this  type  of  design  follows 
the  algorithm:  (1)  design  an  autopilot  which  stabilizes  the 
missile  airframe  and  makes  it  responsive  to  command  inputs, 

(2)  design  a  feedback  controller  which  provides  for  accurate 
tracking  of  the  target  line-o£-slght,  both  stand-alone  and  in 
the  presence  of  missile  maneuvers,  and  (3)  develop  a  guidance 
law  which  guides  the  missile  to  the  target  in  the  desired 
manner.  These  three  tasks -are  traditionally  taken  as  separate 
problems,  and  will  be  treated  as  such  in  this  design. 

Ants  cl  Ian  Design 

The  design  objective  for  the  autopilot  is  to  obtain  a 
stabilized  body  rate  command  system.  This  design  goal  is 
driven  primarily  by  the  guidance  scheme  used  for  proportional 
navigation  described  later  in  this  section.  Autopilot  design 
for  this  airframe  is  complicated  by  two  factors.  First,  the 
missile  is  open-loop  statically  unstable.  This  in  itself  is 
not  a  particularly  difficult  problem,  which  can  usually  be 
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surmounted  by  simple  body  rote  and  normal  acceleration 


feedback.  It  is  complicated  in  this  design  by  nonminimum 
phase  dynamic  response  of  the  missile  normal  acceleration, 
caused  by  tail-mounting  the  fins.  Initially,  a  positive  fin 
deflection,  as  defined  in  Chapter  III,  will  induce  normal 
acceleration  in  the  positive  direction.  However,  the  positive 
fin  deflection  causes  a  negative  angle-of-incidence  rate, 
which  generates  a  negative  normal  force  stronger  than  the 
positive  normal  force  due  to  fin  deflection,  and  a  negative 
normal  force  is  the  final  result,  after  some  time  lag 
dependent  upon  the  specific  airframe  dynamics. 

To  meet  the  design  goal  it  is  necessary  to  design 
feedback  loops  using  sensor  outputs.  The  airframe  dynamics 
and  sensors  are  described  in  Chapter  III.  The  two  outputs 
available  for  stabilizing  the  airframe  are  body  rate  and 
normal  acceleration.  The  design  will  be  accomplished  using 
sequential  loop  closures.  First,  the  body  rate  loop  will  be 
examined  and  closed,  followed  by  the  normal  acceleration  loop. 
Figure  5.1  depicts  the  scheme  used  to  accomplish  the  design 
goal . 


Figure  5.1.  Body  Rate  and  Normal  Acceleration  Feedback 
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The  Transtec  functions  for  these  two  loops  (derived  fr;m 

the  equations  of  motion  in  Chapter  III)  are 

An(S)  77(s+2. 5278) (s-2. 5278) 

<5  (s)  =  (s  +  77)  (s  +  1. 5917)  (s-0. 5717)  (3.13) 

com 

am(s  } _ -211(3+1.1894  ) _ 

<5  (s)  (s  +  77)  (s+1. 5917)  (s-0. 5717)  (3.10) 

Lom 

The  root  loci  for  these  functions  are  shown  in  Figures  5.2  and 
5.3,  respectively.  Note  that  the  scale  of  these  figures 
precludes  inclusion  of  the  locus  for  the  pole  at  s  =  -77. 

This  was  done  to  focus  on  the  action  of  the  critical  poles 
near  the  origin.  Note  also  that  these  loci  are  for  negative 
values  of  the  K's  as  defined  in  Figure  5.1. 


Figure  5.2.  Root  locus  for  Normal  Acceleration  Feedback 
(Negative  Gain) 

To  stabilize  the  missile,  feedback  must  be  used  to  move 
the  system  poles  as  far  as  possible  into  the  left  half-plane 
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while  maintaining  good  damping.  To  make  the  autopilot  a  body 

rate  command  system,  K&  must  be  driven  as  high  as  possible 

m 

(this  is  known  as  increasing  the  loop  gain).  It  is  evident 
from  inspection  of  the  root  loci  that  first  the  gain  of  the 
body  rate  loop  must  be  Increased  to  move  the  unstable  pole 
into  the  left-half  plane.  Also,  it  Is  desirable  to  move  the 
pole  at  -1.5917  until  it  is  left  of  the  A  zero  at  -2.5278  In 
order  to  modify  the  shape  of  the  An  root  locus.  A  negative 
gain  of  magnitude  1  or  greater  will  accomplish  this.  A  body 
rate  gain  of  -2  changes  the  normal  acceleration  loop  transfer 
function  to 

An(3)  77(3+2.5278) (s-2. 5278) 

<5com(s)  '  (s  +  71)  (s  +  5.92)  (s  +  1.03)  (5.1) 

With  the  poles  so  placed,  the  root  locus  for  the  normal 

acceleration  feedback  loop  now  looks  as  in  Figure  5.4. 
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Figure  5.4.  Root  Locus  for  Normal  Acceleration  Feedback, 

with  =  -2.  (Negative  Gain) 
m 

This  locus  is  much  better  than  the  previous  one  from  a  design 
viewpoint.  Evaluation  reveals  that  a  negative  gain  in  this 
loop  is  desired  to  decrease  response  time.  The  leftmost  pole 
shown  in  Figure  5.4  joins  rapidly  with  the  servomotor  pole  at 
s  =  -71  and  this  pair  then  moves  toward  the  right  half-plane 
zero  as  the  gain  is  increased.  Thus,  to  preserve  good 
damping,  a  maximum  gain  of  -0.7  can  be  allowed.  It  is  desired 
to  keep  the  missile  as  stable  as  possible,  so  a  value  of  -0.14 
will  be  used.  This  is  primarily  to  keep  the  gain  and  phase 
margins  [7:430-433]  high,  which  provides  some  robustness  to 
modeling  uncertainties.  The  body  rate  open-loop  transfer 
function  is  given  by 

em(3)  =  _ 422(3  +  1.1894) _  (5#2) 

Q  (s)  s3+67.24s2+77.63s-1.1863 

m 

com 
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The  3ode  phase  and  magnitude  plots  for  this  transfer 
function  are  shown  in  Figure  5.5,  from  which  it  can  be  seen 
that  these  gains  for  the  airframe  control  loop  yield  a  gain 
margin  of  greater  than  40  dB  and  a  phase  margin  of  84  degrees. 

Further  examination  reveals  that  increasing  the  body  rate 
gain  buys  very  little  movement  of  the  rightmost  pole,  while 
bringing  the  other  two  closer  together,  so  that  a  lower 
acceleration  gain  is  allowed  to  maintain  damping  and  preserve 
good  gain  and  phase  margins.  Therefore,  a  body  rate  gain  of 
-2  and  an  acceleration  gain  of  -0.14  provide  about  the  best 
performance  possible. 

The  step  response  due  to  a  body  rate  command  is  shown  in 
Figure  5.6.  This  figure  shows  that  the  design  objective  of  a 
body  rate  command  system  is  met  very  well. 

Tarast  Tracker 

The  design  goals  for  the  tracker  system  are:  (1)  minimize 
target  track  error,  (2)  minimize  response  time,  and  (3) 
minimize  transient  behavior.  For  this  design  effort,  the 
tracker  will  be  considered  detached  from  the  missile.  Once 
again,  successive  loop  closures  will  be  used  to  consider 
feeding  back  both  antenna  rate  and  position.  The  antenna  rate 
loop  will  be  closed  first,  followed  by  the  position  loop.  The 
scheme  used  to  accomplish  this  design  goal  is  depicted  in 
Figure  5.7. 
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Figure  5.5.  0  Open-Loop  Transfer  Function  Bode  Plots 


Figure  5.6.  Airframe  Response  to  Step  Body  Rate  Command 


<t> 


T 


e 
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Figure  5.7.  Tracker  Feedback  Loops 


The  step  response  of  the  'unmodified'  tracker  (with  K^=l 

simply  to  close  the  loop  and  =0)  is  shown  in  Figure  5.8. 

s 

This  response  shows  the  tracker  to  be  slightly  overdamped,  but 

otherwise  a  fairly  good  design  to  begin  with.  As  described  in 

Chapter  III,  all  three  states  of  this  component  are  directly 

available  for  measurement.  However,  only  the  <p  and  ©s 

measurements  are  needed  to  facilitate  the  design,  since  the 

boresight  error  resolver  forms  a  natural  unity  feedback  of  9  , 

s 

as  is  apparent  in  Figure  5.7.  The  transfer  functions  for  the 
loops  of  interest  are  (with  =  0  and  =  1) 


Vs) 

*T(s) 


-1945.8 

s ( s  +  10  4 ) (s+5.26) 


V3> 

w(s ) 


94 

( s  +  10  4 ) 


(5.3) 

(5.4) 


The  root  loci  for  these  transfer  functions  are  shown  in 
Figures  5.9  and  5.10,  respectively.  Looking  at  these  root 
loci,  the  best  course  of  action  is  not  obvious.  It  is 
desirable  to  make  K  .  as  high  as  possible,  to  increase  the  loop 
gain,  but  obviously  this  will  drive  the  poles  rnstable  in 
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Figure  5.8.  Unmodified  Tracker  Response  to  Step  Command 


Figure  5.9. 


Root  Locus  for  K  .  Gain(©  Feedback) 

<P  3 


fcal  s 

Figure  5.10.  Root  Locus  for  9  Feedback 

s 


short  order.  However,  increasing  the  gain  on  will  move 

s 

the  'center  of  gravity'  of  the  loci,  with  favorable  results 

This  can  be  seen  in  the  root  locus  for  K  .  with  Ka  =10.  Th 

4>  e3 

gain  moves  the  servomotor  pole  to  -1044,  and  the  new  K  .  roo 


locus,  shown  in  Figure  5.11,  is  much  more  tractable  from  a 
design  point  of  view. 


Now  K  .  can  be  increased  to  a  value  of  15,  which  provides  high 
loop  gain  and  yet  maintains  good  damping.  The  new  step 
response  of  the  system,  shown  in  Figure  5.12,  displays  this 
behavior.  Examination  of  the  Bode  plot  for  the  open-loop 
transfer  function 

Vs)  ,  mu _  (5  5i 

*T(s)  s ( s  +  2 . 56 ) ( s+1044  )  v  ' 

as  shown  in  Figure  5.13,  reveals  that  the  gain  margin  for  this 

design  is  39  dB,  and  the  phase  margin  is  27  degrees. 

Increasing  Ki  further  would  enhance  the  time  response,  but 
9s 

would  degrade  the  stability  margins,  so  no  further 
modification  will  be  made.  The  resultant  closed-loop  system 
is  shown  in  Figure  5.14. 


5-12 


5-13 


Figure  5.12.  Response  of  Tracker  to  Target  Position  Step  Connand 


Figure  5.14.  Closed-Loop  System  Block  Diagram 
At  this  point,  the  design  would  normally  be  finished. 
Unfortunately,  the  step  response  of  the  coupled  system  (Figure 
5.15)  to  a  body  rate  step  command  shows  that  missile  maneuvers 
are  adversely  affecting  the  tracking  accuracy.  This  is  most 
obvious  in  the  magnitude  of  <p,  the  filtered  boresight  error. 
Also  notice  the  difference  between  and  e_,  which  ideally 
should  be  exactly  opposite  for  perfect  tracking,  since  the 
target  is  not  moving.  This  effect  may  be  minimized  by  feeding 
the  airframe  body  rate  to  the  antenna  rate  command  signal,  <*>. 
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missile  bodjr  rate 
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Figure  5.15.  8yste»  Response  to  Body  Rate  Step  Command 


The  corresponding  block  diagram  for  this  scheme  is  shown  in 

Figure  5.16,  where  is  the  gain  for  the  body  rate  feedback 
to  w,  which  is  to  be  determined. 


Figure  5.16.  Closed  Loop  System  Block  Diagram 


Repeated  step  response  evaluation  reveals  that  a  gain  of  10 
for  Kw  largely  cancels  out  the  transient  effect,  as  shown  in 
Figure  5.17. 


Guidance  Ld*  CoaBiltatlQH 

Traditional  guidance  techniques  for  achieving 
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8*0 


Flqure  5.17.  System  Response  to  Body  Rate  Step  Command 


proportional  navigation  [2:2-14]  indicate  that  the  guidance 
computer  must  command  a  normal  acceleration  of  the  airframe 
proportional  to  the  LOS  rate  of  the  target  relative  to  the 
missile.  The  guidance  lav  is  then  of  the  form 

A  =  (K  ) (LOS  Rate)  (5.6) 

n  Q 

com 

However,  An  is  related  to  the  missile  bearing  angle  y  by 

A  =  y  | V I  (5.7) 

n  11 

and  r  is  given  by 

Y  -  9  -  a  (5.8) 

If  a  is  assumed  small,  then  the  guidance  lav  becomes 

Q  =  (K  ) (LOS  Rate )/ 1 V I  (5.9) 

mg  11 

com  ^ 

Thus,  the  origin  of  the  design  goals  for  the  airframe  and 
tracker  are  made  obvious.  The  constant  of  proportionality  Kg 
for  this  guidance  should  also  be  a  function  of  closing  rate, 
so  the  above  relation  takes  the  form 

9  =  (K  ) |R  | (LOS  Rate)/ | V |  (5.10) 

“com  g 

This  results  in  a  non-linear  problem,  if  both  LOS  rate  and 
closing  rate  vary,  as  they  do  in  a  realistic  encounter. 
Selection  of  the  constant  of  proportionality  is  highly 
dependent  upon  expected  engagement  conditions.  A  lover 
guidance  gain  constant  results  in  pursuit  of  faster  targets, 
and  a  high  one  makes  the  missile  more  responsive  to  target 
maneuvers.  Effectively,  the  guidance  constant  determines  hov 
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long  it  will  take  the  missile  to  correct  its  course  to  the  new 

predicted  intercept  heading.  For  this  missile,  a  value  o£ 

=  2.75  was  chosen  based  on  several  flyout  scenarios  against 
viable  targets.  A  smaller  value  makes  the  missile  fly  a 
pursuit  guidance  on  moderately  fast  (400  ra/s)  targets,  a 
larger  one  causes  it  to  react  too  violently  to  nominal  (6g) 
target  maneuvers.  The  nominal  value  of  (Vj  =  1100  m/s  was 
used  in  the  above  guidance  law  derivation.  Also,  assuming 
that  the  tracker  is  accurately  tracking  the  target,  the  LOS 
rate  is  given  by 

LOS  rate  =  6  +  9  (5.11) 

s  m 

Thus,  the  guidance  law  becomes 

=  2 . 75 | R_ | ( Q  +  £  )/1100 

m T  s  m 

com 

=  0.0025  |Rt|  {0s  +  ej  (5.12) 

During  the  flyout  runs,  the  above  guidance  law  was  found 
to  generate  very  large  initial  accelerations,  due  to  the 
assumption  of  small  a.  Therefore,  an  additional  factor  of  .3a 
was  added  to  the  guidance  lav,  to  limit  the  large  normal 
accelerations.  From  Chapter  III, 

«<s)  *  O (a)  -  r(a)  (3.3) 

B 

which  (using  proper  unit  conversions)  converts  to 

A  ( s ) 

0.(3)  =  ©  (3) - - -  (5.13) 

"  2 

yielding  the  guidance  lav 
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''  l  5  ) 

m 

com 


(5.14) 


=  0.0025  j  r  | \  e  (3)  +  e  (3)  ] 

T  1  m  J 

-  (0.3)6  (s)  -  (0.5) A  ( s ) 
m  n 

Since  the  guidance  problem  is  highly  non-linear,  the  author 
has  decided  to  exclude  this  portion  of  the  problem  from  the 
optimal  design  requirements,  and  will  instead  use  the  same 
guidance  law  on  each  design.  The  above  derived  gains, 
feedbacks,  and  guidance  law  computations  can  be  represented  by 
a  single  feedback  matrix  K  to  be  used  in  the  control  law 


u  =  K  y 
P  P 


(5.15) 


such  that  the  closed-loop  missile  system  system  can  be 
represented  as 


where 


X 

p 

-  Acl 

x  +  B  .  w 
p  cl 

(5.16) 

Yp 

II 

o 

n 

p-* 

x_  +  D  .  w 
p  cl 

(5.17) 

Acl  - 

AP 

pll  p 

(5.18) 

Bcl  = 

Bpl2 

(5.19) 

Ccl  = 

CP 

(5.20) 

Dcl  = 

°P12 

(5.21) 

where  Ap,  Bpn,  Bpl2,  Cp,  DpU,  Dpl2,  xp,  up,  and  yp  are  as 
given  in  Chapter  IV  (Equations  4.13-4.22),  and  K  is  given  by 


K  = 


-1.7+0. 0025*R„ 


-0.29 

10 


0 

-15 


0.0025*Rn 

10 


(5.22) 


This  representation  facilitates  implementation  in  the  digital 
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simulation,  which  is  designed  to  handle  the  state-space 
representation  of  the  missile  system. 
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law 


ZL-  H_2  Ovntr^ller  Synthesis 
The  H2  design  objective  is  to  find  the  feedback  control 

u(s)  =  K  ( s  )  y ( s )  (6.1) 

such  that  the  2-norm  of  the  closed  loop  transfer  function 
matrix  z(s)/w(s)  is  minimized,  where  v(s)  are  the  system 
inputs  (including  commands,  disturbances,  and  noise)  and  z(s) 
is  an  error  signal.  Simply  stated,  it  is  desired  to  minimize 
system  error  response  to  all  signals  entering  the  system. 

The  methodology  for  solving  this  problem  is 
well-established  in  control  theory  literature,  so  a  brief 
overview  will  be  presented  here.  For  a  more  extensive 
treatment,  see  Doyle  et.  al.  [3). 

Theory  lax.  Sialfi.  Space  Solution  oi.  the.  tt2  Design  Pxoblaa 

The  development  presented  in  this  section  was  taken  from 
Doyle  et.  al.  (3J.  The  standard  block  diagram  representation 
of  the  small  gain  problem  is  shown  in  Figure  6.1. 


Figure  6.1.  Standard  Block  Diagram 

In  Figure  6.1,  w  is  the  vector  containing  all  system  inputs,  z 
is  the  vector  containing  error  signals  to  be  controlled,  u  is 
the  control  input  vector,  and  y  is  the  vector  of  measured 
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system  errors.  G  Is  the  open-loop  plant  and  K  is  the 
controller  to  be  designed.  The  resultant  closed-loop  matrix 
of  transfer  functions  from  v  to  z  is  denoted  T 

zw 

2-norm  Definition .  The  objective  of  H 2  design  is  to  find 

K  such  that  T  is  internally  stabilized  and  the  2-norm  of 
zw 

T  ,  denoted  II T  |L,  is  minimized.  The  2-norm  of  T  is 
zw'  11  2Wfl2'  zw 

defined  as 


P 

'i  i. 

lTzvi2  ;= 

1 

In 

“J*00  traCe 

1  T  (jw)*  T  (jw)  1  c 

zw  zw 

J  2 

(6.2 

• 

-oo 

i.  j 

J 

where  the  superscript  *  denotes  the  complex  conjugate 
transpose . 

This  concept  can  be  visualized  easily  in  the  single-input 
single-output  (SISO)  case,  where  the  2-norm  is  proportional  to 
the  area  under  the  Bode  plot  of  the  closed-loop  transfer 
function . 

Visualization  in  the  multi-input  multi-output  (MIMO)  case 
is  more  difficult,  but  the  principle  is  the  same.  In  the  MIMO 
case,  the  2-norm  is  proportional  to  the  sum  of  the  areas  under 
all  of  the  singular  value  plots.  Singular  values  are  used  to 
extend  the  concept  of  bode  magnitude  plots  to  the  MIMO  case. 
For  an  extensive  treatment  of  the  properties  of  the  singular 
value  plot,  see  AFWAL  TR-85-3102  [81. 

Problem  Requirements  and  Scaling .  In  state  space  form,  G 
can  be  represented  as 


x=Ax+B1w+B2u 

(6.3) 

z  =  Cl  x  *  Du  w  ♦  D12  u 

(6.4) 
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Y  :  C2  X  f  D21  W  f  D22  U 

To  solve  the  problem  via  existing  theory  133,  the  following 
conditions  must  be  satisfied; 


1) 

DU 

=  0 

(6.6) 

2) 

°22 

=  0 

(6.7) 

3) 

(a ) 

D12T  d12  *  1 

(6.8) 

(b) 

°21  D21  =  1 

(6.9) 

4) 

(A, 

B^)  stabilizable  & 

(C1#A) 

detectable 

5) 

(A, 

)  stabilizable  & 

(c2,a) 

detectable 

Condition  1  requires  that  the  vector  of  error  signals  z 
receives  no  direct  feedforward  from  the  disturbances  w.  If 
this  requirement  is  violated,  ST  I.  is  equal  to  oo  for  any 

"  ZW  "  L 

K ( s )  developed. 

Condition  2  requires  that  the  vector  of  measured  errors 
receive  no  direct  feedforward  of  the  control  inputs  u.  This 
condition  can  actually  be  relaxed  if  required  [51  &  [111,  but 
this  is  not  necessary  for  this  problem. 

Condition  3a  requires  that  z  receives  direct  feedforward 
of  each  element  of  u,  and  Condition  3b  requires  that  each 
element  of  y  receives  direct  feedforward  of  some  element  of  w. 

In  summary,  each  measured  output  must  be  directly 
contaminated  by  some  disturbance  but  not  by  any  control  inputs 
(as  will  be  the  case  in  any  well-modeled  system)  and  the  error 
signal  z  must  be  affected  by  each  control  input  directly 
(every  control  input  must  be  controlled)  but  not  directly  by 
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any  of  Cue  J  btar  bance  inputs  {penalizing  tiie  J  iotut  bance 

Inputs  directly  would  yield  an  infinite  2-norm  since  they 
can't  be  affected). 

Conditions  3a  and  3b  also  require  that  the  input  u  and 
the  output  y  be  scaled.  This  scaling  is  accomplished  via  a 
Cholesky  decomposition  [3],  which  yields  unique  upper 
triangular  matrices  that  satisfy 


where  the  superscript  -1  Indicates  inverse  and  the  superscript 
-T  indicates  transpose  of  the  inverse.  The  scaling  is 
appended  to  the  original  plant  G  as  shown  in  Figure  6.2,  where 
the  superscript  ~  denotes  scaled  values. 


Figure  6,2.  Block  Diagram  of  Plant  with  Scaling 


The  scaling  will  be  considered  part  of  the  plant  in 
determining  the  controller,  but  will  be  considered  part  of  the 
controller  after  the  solution  Is  obtained.  Appending  the 
scaling  to  the  plant  G  results  in  a  modification  to  the 
original  plant  representation.  These  modifications  can  be 
represented  as  follows,  where  the  superscript  ~  denotes  the 
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scaled  matrices,  which  are  given  by 


B.  =  B-  S  1 

(6.12) 

2  2  u 

C2  =  Sy  C2 

(6.13) 

~  -1 

D12  =  °12  Su 

(6.14) 

°21  =  Sy  °21 

(6.15) 

The  problem  is  now  in  the  form  required  for  application  of  the 
existing  methods  for  derivation  of  the  controller,  as 
described  in  Doyle  et.  al.  131. 

Solution .  This  method  requires  obtaining  the  real. 


unique. 

symmetric  solutions  X2  and  Y2 

to 

the 

two  uncoupled 

Riccati 

equations 

F1T  X2  *  X2 

F1  -  X2  G1  X2  + 

H1 

=  0 

(6.16) 

f2t  y2  ♦  y2 

f2  -  y2  g2  y2  ♦ 

H2 

=  0 

(6.17) 

where 


F,  = 


A  B2  D12  C1 


G,  = 


H.  = 


B2  B2 

-  T  - 
C!  C1 


C,  = 


(  1  -  D12  d12*)  C, 


F2  "  A  '  B1  °21  C2 
g2  -  a2T  c2 

A*  A 

H2  =  B1  B1 


B,  » 


B1  (  1  ~  D21  °21> 


(6.18) 

(6.19) 

(6.20) 
(6.21) 
(6.22) 

(6.23) 

(6.24) 

(6.25) 


The  optimal  controller  (denoted  by  the  subscript  cp)  is  then 
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given  by 


B.  K 


K 


£ 


0 


where 


(6.26) 

(6.27) 

(6.28) 
(6.29) 


Kc  =  B2T  X2  t  5UT  (6.30) 

K£  =  Y2  J2T  V  Bl  S21T  (6.31) 

Problem  Set-up 

The  problem  now  standing  in  the  way  of  deriving  a 
controller  is  the  conversion  o£  the  missile  model  into  one 
conforming  to  the  standard  problem. 

The  design  objective  is  to  derive  a  controller  such  that 
the  antenna  will  tr  ;ck  the  target  and  the  missile  can  be 
steered.  This  is  a  tracking  problem  and  Francis  [4]  presents 
a  methodology  to  convert  the  tracking  problem  to  the  small 
gain  problem.  That  methodology  will  be  employed  here. 

The  Tracking  Problem.  To  convert  the  tracking  problem  to 
the  small  gain  problem,  it  is  necessary  to  convert  the  desire 
to  accurately  track  the  target  and  guide  the  missile  into  a 
requirement  to  minimize  some  set  of  values.  It  is  desired  to 
minimize:  (1)  error  between  commanded  and  actual  missile  body 
rate,  (2)  transients  in  normal  forces,  (3)  error  between 
antenna  boresight  and  target  1 lne-o£-sight,  (4)  transients  in 
antenna  rates,  and  (5)  control  inputs.  Expressing  these 
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desires  in  terms  of  the  tracking  problem  presented  in  [41 

yields 


where 


z 


r-v 

.  pu. 


Y  -  [r-v] 


(6.32) 

(6.33) 


p  =  diagonal  matrix  of  control  weighting  factors 
r  =  reference  signal  desired  to  track 

v  =  controlled  variables,  which  track  reference  signal 
u  =  control  inputs 

For  the  missile  system  these  become 


v  = 


(missile  body  rate) 
(normal  acceleration) 
(boresight  error) 
(antenna  rate) 
(antenna  position) 


(6.34) 


r 


(missile  body  rate  command) 

(desire  no  acceleration  transient) 
(desire  no  boresight  error  transient) 
(desire  no  antenna  rate  transient) 
(antenna  position  command) 


(6.35) 


u  = 


(fin  deflection  command) 
(tracker  command) 


(6.36) 


In  block  diagram  form,  the  system  looks  as  shown  in 
Figure  6.3. 


ll 


Q_J- 


:: _ L  +L 

G - 


->  v 


Figure  6.3.  Converted  System  Block  Diagram 

Conversion  to  Small  Gain  Problem.  It  is  desired  to  express 
the  above  system  in  the  form 

x=Ax+Blw+B2u  (6.3) 


z  =  Cl  x  +■  Dll  v  +  D12  u 
y  =  C2  x  +  D21  v  +  D22  u 


(6.4) 

(6.5) 


where  K  is  considered  'open'  and  v  Is  a  vector  containing  all 
external  inputs.  The  external  inputs  will  be  arranged  as 

r 

w  =  I  *T  (6.37) 

? 

T? 

This  yields,  in  terms  of  the  system  representation  given  in 
Chapter  IV, 


A  ■  AP 

(6.38) 

81  -  1  0  Bpl2  Bp2  0  1 

(6.39) 

82  ■  Bpll 

(6.40) 

Cl  * 

ttf 

(6.41) 

C2  = 

[-S] 

(6.42) 

Dll  = 

_I  -Dpl2  "Dp2  -1  ‘ 
0  0  0  0 

(6.43) 
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D 1  2 


6  .  i  4  i 


Figure  6.4.  System  with  Weights 


In  this  system,  each  of  the  weighting  functions  is 
frequency  dependent.  When  expressed  in  state  space  form,  each 
weighting  function  generates  additional  states,  one  for  each 
pole  present  in  the  transfer  function  describing  it. 

Therefore,  each  weighting  function  has  its  own  A,  B,  C,  and  D 
matrices.  Denoting  the  respective  weighting  function  matrices 
with  the  subscripts  associated  with  the  W's  in  Figure  6.4 
yields  the  new  set  of  matrices  describing  the  total  system 
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(with  K  open  as  be£oce) 


B1  = 

B2  = 

Cl  = 

C2  = 

Dll  = 

D12  = 

D21  = 


0 

0 

0 

-B  C 
z  p 

'  0 

B1 

0 

0 

B  D. 

Z  1 


0 

A1 

0 

0 

B  C, 
z  1 


[Bpl2  Bp21C3 
0 

A_ 


-B  CD  . .  D  » ]C. 
z  pl2  p2  3 


0 

0 

0 

A3 

■BzC2 


(Bpl2  Bp21D3 
0 

B_ 


-BzCDpl2  Dp21D3 


0 

0 

0 

B  3 

-BzD2 


-DzDpll 

P 


[D1  ~°2  -lDpl2  DP21D3] 


0 

0 

0 

0 

A 


(6.48) 


V 

0 

0 

0 

BzDpll 


DzCp  DzCl  *DzC2  -Dz[DP12  DP2]C3  Cz 


[  Cp  C1  C2  CDpl2  Dp21C3  °] 


[DzDl  -DzD2  'Dz1Dp12  Dp21D3 


(6.49) 


(6.50) 


(6.51) 


(6.52) 


(6.53) 


(6.54) 


(6.55) 
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D  2  2 


(6.56) 


Now,  Condition  1  translates  to  D  =0.  Condition  3a  is  met 
(using  scaling)  as  long  as  all  the  control  input  weights  are 
non-zero,  and  3b  is  met  as  long  as  care  is  taken  that  and 

D2  are  chosen  such  that  each  signal  in  y  receives  some  direct 
feed-forward  from  at  least  one  of  the  disturbances  (#T,  ?,  or 
t?)  • 


In  addition  to  satisfying  Condition  1,  the  weighting 
functions  are  a  powerful  design  tool.  They  enable  the 
designer  to  shape  the  system  response.  This  is  accomplished 
by:  (1)  designing  the  filter  on  the  incoming  disturbances  such 
that  the  Bode  plot  of  the  filter  looks  like  the  Fourier 
Transform  of  the  incoming  signal;  and  (2)  designing  the  filter 
on  the  error  signal  such  that  the  desired  portion  of  the 
spectrum  of  the  error  is  penalized. 

FLiLeElag  the  Disturbance  Input.  The  purpose  of 
designing  the  disturbance  input  filters  is  to  properly  reflect 
the  frequency  content  of  the  incoming  signals.  The  automated 
state  space  solutions  to  the  problem  are  constructed  to 
minimize  the  2-norm  of  Tzw«  At  this  point,  it  is  necessary  to 
remember  that  the  transfer  function  of  a  system  is  equivalent 
to  the  impulse  response  of  the  system.  Therefore,  minimizing 
the  2-norm  of  Tzw  is  equivalent  to  minimizing  the  error 
response  to  an  impulse  input.  In  reality,  the  system  will  not 
see  Impulse  inputs,  so  minimizing  the  2-norm  of  T  as  is  puts 
an  unwanted  requirement  on  the  design.  To  express  the  true 
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design  goals  ,  it  is  necessary  to  convert  T  to  the  true 

zv 

system  response  desired. 

As  an  example,  consider  the  case  of  a  command  w  which 
enters  the  system  as  a  step  input.  The  Fourier  Transform  of  a 
step  function  is  Plotted  as  magnitude  (in  dB)  versus 

frequency,  the  transform  looks  as  shown  in  Figure  6.5. 
Obviously,  a  step  function  has  very  low  magnitude  at  high 
frequency  and  large  magnitude  at  low  frequency.  Now  consider 
a  SISO  plant  P  to  which  this  command  is  applied  as  shown  in 
Figure  6.6,  where  y  is  the  response  and  z  is  the  error  to  be 
minimized . 


w 


z 


Figure  6.6.  SISO  Plant  Block  Diagram 


Two  cases  will  be  examined,  for  comparison  purposes. 
Representing  these  as  PI  and  P2,  let 

P1(J">  *  -jsn-  (6-51) 

P2(J“)  *  -JSTTJ-  (6-58> 

which  yield  the  time  responses  shown  in  Figures  6.7  and  6.8, 
respectively.  These  systems  have  time  constants  of  1  and  10, 
respectively,  and  both  have  zero  steady-state  error. 
Accordingly,  the  quality  of  these  systems  should  be  reflected 
in  the  size  of  the  2-norm  of  T  .  For  plant  PI, 

Z  V 


T  ( jw) 
zw 


z  (  jw) 
w(  jw) 


l-Pl(ju) 


jw 
j  w+1 


(6.59) 


6-12 


1.2 
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Figure  6.7.  Step  Response  of  PI 
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Taking  the  2-norm, 


zw 11 2 


1 

p+00 

trace  | 

jw  * 

jw 

•> 

_  1  A,  v 

2n 

J 

-oo 

jw+1 

jco+1 

—  I  Cl  OJ 

1 

p+OO 

trace 

J  dw 

J 

1 

2 

2* 

w 

J 

-oo 

L  l+CO2 

- 

1 

L  - 

-  tan  1(w)J 

00 

1 

2 

2  rr 

L 

-oo  . 

=  00 


Similarly,  for  P2, 


T2v(j“'  =  -Invj-  =  1-p2<^'  ■ 


Taking  the  2-norm, 
ITzJ2 


1  r*®  trace  f 

jco  * 

jw 

|  dco 

2ft  J  tEaCe  1 

jw+10 

j  u>+ 1 0 

¥ 

(6.60) 


(6.61) 


1 

2 


=  yn-  J^00  trace  [  - 1  do> 

2n  -oo  «-  100+co  J 


1 

2 


(-M 


w  -  10  tan 


1<W10,]|1)T 


=  oo  (6.62) 

Clearly,  this  computation  does  not  yield  a  very  good  measure 
of  actual  system  performance  for  comparison  purposes.  Again, 
this  is  because  the  T  examined  here  is  the  error  response  to 

2  V 

an  impulse  input.  The  effect  of  adding  a  filter  with 
frequency  characteristics  matching  that  of  a  step  function 
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will  aow  be  examined. 

The  Fourier  Transform  of  a  step  function  is  l/jco.  Adding 
a  filter  with  this  frequency  characteristic  yields  a  new 
system  as  shown  in  Figure  6.9. 


figure  6.9.  SISO  Plant  with  Filter 


Now  T  for  PI  looks  like 
zw 


zw 


iT,J 


zw  11 2 


_ 1 

jw 

1 

1- j  to 

1  + jtO 

1  + j  to 

l+to2 

-norm  . 

as  before. 

* 

1 

J400  trace  | 

■  1-jw  * 

1- jto 

In 

•  1+W2 

l  +  u>2 

1 

J400  trace 
-00 

r  i 

]  d“) 

In 

* 

L  i+«2 

(6.63) 


d» 


1 
2 


1 

I 


■  (-h-  r 


■( 


1 

2 


2n 


rr 


(6.64) 


For  P2, 


Tzw(jw)  " 


1  j  0) 

j  to  10  + jw 


1 

10  + jw 


1 0  -  J  q> 
100+w 


(6.65) 


6-17 


Computing  the  2-norm 


2 


1 

2  n 


j+QD 

-00 


trace 


[ 


1 0  -  j  co  * 
100+co2 


10-  jco  "J 

100+co2  J 


dco 


L 

2 


- 

1 

J+°°  trace  f 

i 

|  dw 

2n 

,  ^  2 

s 

-oo  L 

100+w  - 

* 

1 

r  i 

tan  1(w  /  10)J 

00 

2n 

» 

L 

-00  , 

* 

1 

n  /  10 

1  ? 

2  TT 

b 

J 

m 

10 


(6.66) 


Clearly,  this  is  a  much  better  reflection  of  the  true 
performance  of  the  systems  in  response  to  a  3tep  input.  The 
key  is  to  tailor  the  disturbance  filter  to  look  like  the 
expected  disturbance  (or  command)  input.  To  get  an  even 
better  idea  of  how  close  the  system  is  responding  to  the 
desired,  it  is  necessary  to  examine  filtering  of  the  error 
signal . 

Filtering  Si.  Ltlfi.  Error  Signal.  Suppose,  using  the 
plants  above,  that  it  is  desired  to  penalize  only  the  steady 
state  error  of  the  system.  Since  the  above  system  has  no 
steady  3tate  error,  proper  filtering  of  the  error  should 
reflect  this.  If  it  Is  desired  to  penalize  steady  state  error 
only,  then  the  gain  of  the  filter  should  be  high  on  the  low 
frequency  end  and  low  on  the  high  frequency  end.  A  filter 
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which  reflects  this  desine  is  simply  1/ (  j 1  )  .  Appending  this 

filter  to  the  error  output  z  yLelds  the  system  3hown  in  Figure 
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Figure  6.10.  SISO  Plant  with  Filters  On  Input  and  Output 


The  transfer  function  T  for  PI  is  then  given  by 

zw 


.  _  _ 1_  jco  1 

Ztfi:,U>  j<*>  1  +  jco  jcO+1 


and  the  2-norm  of  T  becomes 

zw 


( l-c*>  )-2joo 
2  2 

(!+«■  ) 


(6.67) 


|T 


(  1-oj2  )  -2  jw  *  (l-w2)-2j« 


22 

(1+w  > 


h  =  trace  [ 

■(-fc-C  [-mV]-) 


TT 

(1+w  P 


H 


1 

2 


1 

2 


■(-fed 

■(-fed 


CO 


2 ( w  +1 ) 


(  0  + 


* 


)  -  (  0 


1 

2 


1 

2 


For  P2,  the  transfer  function  becomes 


T  t  -  1  1  (  10  -<1>  )  -11  jg> 

zw  ~jzr  nti»  uoo*«2)(it<w5) 


(6.68) 


(6.69) 
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and  the  2-norm  becomes 


!Tzv»2 


1 

2 


1  j-  r  - ^ —  d«  1 1 

2  —oo  L  (  100+w  )  ( 1+'A>2 )  j  J 

-57 -  [  4r-(-n-tanl“/10)tta"‘1|“>]]La, 


•  '  4r  [  M'S-  "  H  ]  ) 


1 

2 


1 

2 


■  hV) 


1 

2 


=  0.067 


(6.70) 


Now  the  true  nature  of  this  process  is  revealed,  as  the  norm 
of  p2  is  very  much  smaller  than  that  of  PI,  due  to  the  fact 
that  the  error  resultant  from  P2  has  much  less  low  frequency 
content  (faster  response)  than  that  for  PI.  The  principles  of 
behavior  demonstrated  in  this  SISO  example  extend  directly  to 
the  MIMO  problem,  with  the  added  concept  of  relative  weights 
between  error  signals. 

Relative  Weighting  ai.  GLaaponenta  ai  the.  Error 
Signal .  In  composing  the  problem,  the  ability  to  modlfv  the 
relative  Importance  between  the  elements  of  the  error  signal  z 
allows  the  designer  to  tune  the  response  to  his  desires  when 
some  of  the  signals  compete  with  each  other.  The  most  obvious 
competition  takes  place  between  the  control  signal  u  and  the 
error  signal  z,  as  decreasing  the  penalty  on  u  allows  more 
control  to  be  used  to  decrease  z. 

In  a  MIMO  system,  the  same  competition  can  occur  between 
elements  of  the  error  signal  itself.  An  obvious  case  of  this 


6-20 


occurs  if  the  desire  is  to  penalize  some  tracking  error  and  at 
the  same  time  penalize  the  rate  of  the  subject  response.  it 
is  not  possible  to  maintain  a  zero  rate  and  still  respond  to 
the  command.  The  determination  of  the  relative  penalties  is 
almost  always  an  iterative  process.  If  the  response  is  too 
slow,  the  position  error  is  penalized  more  heavily.  Likewise, 
if  the  response  is  too  lightly  damped,  the  rate  is  penalized 
more  heavily. 

In  this  context,  the  flexibility  of  the  weighting  filters 
shows  its  true  potential.  With  the  addition  of  the  frequency 
dependent  weights,  not  only  can  the  designer  set  the  relative 
penalties  on  the  different  errors  overall,  he  can  also  set 
relative  penalties  between  errors  that  vary  over  the  entire 
spectrum.  In  the  case  of  the  position  versus  rate  errors, 
this  flexibility  allows  the  designer  to  penalize  the  position 
more  at  low  to  intermediate  frequency,  so  that  the  steady 
state  error  and  rise  time  are  reduced,  and  simultaneously 
penalize  the  rate  in  some  frequency  range  where  the  system  is 
exhibiting  poor  damping.  This  was  the  approach  used  in 
developing  the  weights  on  the  missile  system. 

The  weighting  function  W  can  be  represented  as 

Z 

*z  -  Pz  *z  l6-111 

where  p  is  a  diagonal  matrix  containing  the  relative  weights 
z 

between  the  components  of  the  error  signal  and  Wz  is  a 
normalized  frequency  dependent  weighting  function  which 
achieves  unity  gain  at  its  maximum  value.  This  is  done  to 
separate  the  magnitude  of  the  filtered  signal  from  the  shape 
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of  the  filter,  facilitating  ease  of  design.  By  making  such 
a  separation.  It  is  possible  to  automate  the  calculation  of 
W  whenever  the  weighting  factors  p  or  the  shaping  filters  W 

z  Z 

are  modified. 

All  of  the  concepts  needed  to  begin  the  design  have  been 
developed.  It  quickly  becomes  evident,  however,  that  with  the 
many  design  variables  available  (shape  of  input  and  output 
weighting  filters,  relative  weights  between  error  signals  and 
control  inputs,  etc.)  that  the  design  must  be  carried  out  in 
some  systematic  fashion,  or  else  the  designer  can  be 
overwhelmed  by  the  many  options  available.  An  algorithm  was 
therefore  developed  to  simplify  this  process. 

Design  Procedure .  The  first  step  in  the  design  process 
is  to  break  the  tasking  into  manageable  objectives.  The 
design  task  was  therefore  broken  into  three  portions:  (1) 
target  tracking,  (2)  missile  command  following,  and  (3) 
disturbance  filtering. 

Target  Tracking .  First,  all  input  weighting  filters 
were  set  to  identity.  Next,  all  weights  on  the  body  rate  and 
normal  acceleration  error  were  set  to  zero,  to  focus  on 
tracking  the  target  only,  and  the  weighting  on  6  was  reduced 
to  a  small  magnitude. 

At  this  point,  it  was  found  to  be  necessary  to  raise  the 
penalty  on  fin  deflection,  <5,  to  the  same  magnitude  as  the 
weight  on  the  tracker  command  signal,  o>,  because  the  missile 
was  attempting  to  maneuver  the  whole  missile  body  in  order  to 
align  the  antenna  with  the  target  due  to  the  low  penalty  on  <5. 
This  is  shown  in  Figure  6.11,  where  the  input  is  a  ramp  of  one 
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Figure  6.11.  System  Response  to  Ramp  and  Step  Target 


decree  pet  second  in  the  target  position  input. 

Successive  Iterations  revealed  that  the  weighting  mix 
that  achieves  the  best  target  tracking  available  with  this 
configuration  was 


p 


z 


0  0  0  0  0 

0  0  0  0  0 

0  0  0.1  0  0 

00  0  0.1  0 
0000  10000 


p 


0.1 

0 


0 

0.1 


(6.72) 


(6.73) 


This  yields  the  response  shown  in  Figure  6.12  for  a  unit  ramp 
antenna  position  command  input  followed  by  a  negative  unity 
step.  From  this  figure,  it  can  be  seen  that  the  transient 
response  of  the  tracker  is  underdamped,  and  a  steady  state 
error  of  25%  is  present  on  the  ramp  input.  To  overcome  this 
error  and  transient,  it  is  necessary  to  Introduce  the 
frequency-dependent  weighting  on  both  the  input  command  w1  and 
the  output  z. 

To  reflect  the  nature  of  the  expected  input  commands  to 
antenna  position,  the  filter  is  modified  (from  the  flat 
spectrum  unity  gain  used  at  the  start)  such  that  the  weighting 
on  antenna  position  commands  looks  as  shown  in  Figure  6.13, 
simulating  a  step  command  over  the  frequency  range  of 
interest.  Also,  the  weighting  on  the  antenna  position  error 
signal  13  modified  to  make  it  a  low  pass  filter  with  a 
bandwidth  of  40  radians  per  second,  as  shown  in  Figure  6.14. 
The  weighting  filter  on  the  antenna  rate  was  modified  to  look 
like  a  bandpass  filter,  as  shown  in  Figure  6.15.  In  addition. 
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figure  6.12. 


lgure  6.14.  Antenna  Position  Error  Weighting  Filter 


figure  6.15.  Antenna  Rate  Error  Weighting 


the  weighting  matrix  on  the  error  signals  was  modified  to 


P 


z 


0  0  0  0  0 

0  0  0  0  0 

0  0  0.1  0  0 

0  0  0  10  0 

0  0  0  0  1000 


(6.74) 


The  response  of  the  resultant  controller  is  shown  in  Figure 
6.16.  This  response  is  satisfactory.  With  the  tracker 
designed,  the  missile  control  system  is  the  next  subject  of 
interest . 

Missile  Control  System.  The  above  process  was 
applied  to  the  missile  command  signals  and  errors.  This 
process  required  modifying  the  body  rate  command  filter  to 
look  like  a  step  function  over  the  frequency  range  of 
interest,  as  shown  in  Figure  6.17.  The  body  rate  error  filter 
was  modified  to  look  like  a  lowpass  filter  to  penalize 
steady-state  error  more  heavily  than  high-frequency  error,  as 
shown  in  Figure  6.18.  In  the  same  manner  that  antenna  rate 
error  was  bandpassed  to  reduce  the  antenna  transient,  the 
normal  acceleration  error  filter  was  converted  to  a  bandpass 
to  achieve  the  same  effect,  as  shown  in  Figure  6.19.  The 
weighting  matrix  on  the  error  signals  was  modified  to 


(6.75) 


These  modifications  yielded  the  response  shown  in  Figure  6.20 
for  a  missile  body  rate  command  of  1  degree  per  second.  This 
response  is  acceptable.  Now  that  the  design  responds 
correctly,  it  only  remains  to  shape  the  disturbance  inputs. 
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Figure  6.16.  Tracker  Response  to  Raap  and  step  Target  inputs 
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Igure  6.16. 
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Figure  6.19.  Normal  Acceleration  E*  :or  Weighting  Filter 
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Figure  6.20.  Airfraae  Response  to  Ranp  and  step  Target  Inputs 


disturbance  weighting  filters  developed  for  the  missile  syste 
are  described  in  Table  6.3. 


Table  6.3.  Disturbance  Input  Weighting  Filters 

Disturbance 

Filter 

Descr iption 

Body  Rate  Command 

-20  dB/decade  above  0.1 
rad/sec,  low  frequency  gain 

20  dB,  to  simulate  a  step 
command  at  high  frequency. 

Normal  Acceleration 
Command 

Flat  at  0  dB  over  entire 
spectrum.  Normal  acceleration 
is  used  for  damping  only, 
no  command . 

Boresight  Error 
Command 

Same  as  Normal  Acceleration 
Command . 

Antenna  Rate 

Command 

Same  as  Normal  Acceleration 
Command . 

Antenna  Position 
Command 

-20  dB/decade  above  0.01 
rad/sec.  Low  frequency  gain 

58  dB, simulate  step  command. 

Accelerometer 

Noise 

High  pass  above  20  rad/sec, 

-40  dB  at  low  frequency. 

Rate  Gyro 

Noise 

High  pass  above  10  rad/sec, 

-40  dB  at  low  frequency. 

Potentiometer 

Noise 

High  pass  above  40  rad/sec, 

-40  dB  at  low  frequency. 

Turbulence 

High  pass  above  1  rad/sec, 

-40  dB  at  low  frequency. 

Glint 

High  pass  above  50  rad/sec, 

-40  dB  at  low  frequency. 

Power  Supply 

Noise 

High  pass  above  400  rad/sec, 

-40  dB  at  low  frequency. 

See  Appendix  A  for  the  respective  transfer  functions.  Bode 


P  i '  j  -  w  ,  ^  *  d  o  c  a  -  o  p  a  c  e  .  Inclusion  of  £ncs^ 

weighting  filters  required  further  modification  to  the  command 
and  error  filters.  These  modifications  are  reflected  in 
Table  6.4. 


Table  6.4.  Error  Signal  Weighting  Filters 

Error  Signal 

Filter  Description 

Body  Rate 

Lowpass  below  1  rad/sec, 
dropping  off  at  -20  dB/decade 
thereafter . 

Normal  Acceleration 

0  dB  over  the  range  30  to  190 
rad/sec . 

Boresight  Error 

Lowpass  below  0.01  rad/sec, 

-20  dB/decade  thereafter 

Antenna  Rate 

Error 

0  dB  over  the  range  10  to 

200  rad/sec. 

Antenna  Position 

Error 

0  dB  out  to  40  rad/sec 
dropping  off  at  -20  dB/decade 
thereafter . 

Again,  see  Appendix  A  for  a  detailed  description  of  the 
filters  described  in  Table  6.4. 

The  set  of  weights  used  in  the  final  design  is  given  by 


P 


z 


400  0  00 

0  0.1  0  0 

0  0  0.1  0 

0  0  0  10 

0  0  0  0 


0 

0 

0 

0 

1000 


(6.76) 


This  set  of  weights  and  filters  yielded  a  compensator  which 
responds  as  shown  in  Figure  6.21. 


ClQSSd-LQQP  System  Representation. 

Once  the  controller  is  obtained,  it  can  be  placed  into 
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an  [BA 
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Figure  6.21.  Systen  Response  to  Ra«p  and  Step  Target  Input 


the  original  system,  so  that  the  ciosed-ioup  ayatem  n 

represented  as  shewn  in  Figure  6.22,  where  [A  B  C  1 

cp  cp  cp 

correspond  to  the  designed  compensator  and  (ABC] 

P  P  P 

correspond  to  the  original  system.  Note  that  the  D  matrices 
for  the  plant  and  the  compensator  have  been  omitted  for 
clarity,  since  they  are  all  identically  zero.  is  the 

matrix  which  derives  the  error  signals  from  the  sensor  values. 
Creation  of  requires  implementing  the  proportional 
navigation  laws  described  in  Chapter  V. 


Figure  6.22.  Closed  Loop  System  Block  Diagram 
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(6.73) 

(6.79) 

(6.80) 

(6.81) 

(6.82) 

(6.83) 

(6.84) 

(6.85) 

(6.86) 


Thus,  the  closed  loop  system  is  derived  directly  from  the  open 
loop  description  and  the  designed  compensator. 


Compensator  Order  Reduction 

Due  to  the  nature  of  the  Riccati  equations  (6.16  &  6.17), 
the  compensator  obtained  in  the  above  procedure  is  of  the  same 
order  as  the  augmented  plant;  23  states.  When  added  to  the 
plant,  the  resulting  closed-loop  system  has  a  total  of  30 
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3  faces.  Thi.5  oystem  .3  tar  too  to  simulate  m  an/ 

efficient  manner  In  PC-MATLAB.  Due  to  array  size  limitations 
and  time  step  requirements,  it  was  desired  to  reduce  the 
compensator  to  no  more  than  thirteen  states. 

The  first  attempt  at  reduction  was  to  perform  a  minimal 
realizati  ,»  of  the  system,  to  eliminate  any  pole-zero 
cancellations .  This  was  accomplished  via  the  PC-MATLAB 
function  MINREAL,  which  allows  the  user  to  specify  a  tolerance 
for  pole-zero  cancellations .  For  a  tolerance  of  0.01,  two 
states  were  removed.  For  a  tolerance  of  0.1,  three  states 
were  removed,  but  a  detrimental  effect  on  system  accuracy  was 
observed.  Accuracy  was  evaluated  by  feeding  a  random  signal, 
a  step  command,  and  a  ramp  command  to  both  the  full  order 
compensator  and  the  reduced  order  compensator,  then  examining 
the  output  control  signals.  At  this  point,  it  became  obvious 
that  the  minimal  realization  would  not  accomplish  sufficient 
order  reduction  without  serious  loss  of  accuracy,  so  model 
order  reduction  via  Hankel  singular  value  elimination  [10] 
vas  attempted. 

Using  this  method,  the  compensator  was  reduced  to  11 
states  without  seriously  compromising  efficiency.  This  was 
sufficient  to  make  simulation  possible.  Both  the  full  order 
compensator  and  the  reduced  order  compensator  are  described  in 
Appendix  C. 
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«’  I  l.t.  ii  C  v  M 1 1  w  1 1  £  I  3  yntiicw 

CO 

The  H  design  objective  is  to  find  the  feedback 
00 

compensator  K(s)  such  that  the  oo-norm  of  the  closed  loop 
transfer  function  matrix  z(s)/w(s)  is  minimized,  where  v(s) 
and  z(s)  are  as  described  in  Chapter  VI.  The  methodology  for 
solving  this  problem  is  a  relatively  new  development,  so  a 
more  complete  treatment  will  be  given  than  in  Chapter  VI.  For 
a  complete  treatment,  see  Doyle,  et.  al.  [3]. 

T,h.e.ai.y  tax.  state,  ajaate.  solution  ai.  the.  teals n  Problem 

The  development  presented  in  this  section  was  taken  from 
[31.  The  standard  block  diagram  representation  of  the  small 
gain  problem  is  shown  in  Figure  7.1.  Notice  that  the  same 
basic  setup  is  used  as  for  the  H2  theory. 


Figure  7.1.  Standard  Block  Diagram 

As  in  Chapter  VI,  w  is  the  vector  containing  all  system 

inputs,  z  Is  a  vector  of  error  signals  to  be  controlled,  u  is 

the  control  input  vector,  and  y  is  the  vector  of  measured 

system  errors.  The  resultant  closed-loop  matrix  of  transfer 

functions  from  w  to  z  is  denoted  T  . 

zv 

oo-norm  Definition.  The  objective  of  design  is  to  find 

K  such  that  T  is  internally  stabilized  and  the  oo-norm  of 
zw  1 
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x  O 


minimi  zed  . 


The  'O-norm  of  T 


T  ,  denoted  ;jT  ;]  , 

ZW  11  Z  W  11  K 

defined  as 

||T  II  :  =  sup  o  T  ( jw) 

H  zw11®  max  j_  zw  J 

where 


zw 


i  5 


(7.1) 


a  :=  maximum  singular  value  (7.2) 

max 

The  above  is  read,  "the  no-norm  of  T  is  defined  as  the 

zw 

supremum  over  all  frequencies  of  the  maximum  singular  value  of 
T  This  concept  can  again  be  visualized  easily  in  the  SISO 

case,  where  the  oo-norm  is  equal  to  the  maximum  value  of  the 
Bode  plot  of  the  closed-loop  transfer  function.  In  the  MIMO 
case,  substitute  "maximum  singular  value  plot"  for  "Bode  plot" 
in  the  above.  Thus,  the  ao-norm  is  a  measure  of  the  worst  case 
amplification  (in  the  frequency  domain)  of  the  output  z  in 
response  to  the  input  w,  in  contrast  to  the  2-norm  being  a 
measure  of  average  amplification  over  the  entire  frequency 
spectrum.  For  this  reason,  controller  designs  based  on 
optimizing  the  ao-norm  return  controllers  which  make  T  (jw) 
have  a  maximum  singular  value  plot  which  is  fiat,  i.e.  all 
frequencies  are  amplified  at  the  same  level. 

Like  the  design  method,  solution  of  the  H  problem 
rests  on  solving  two  uncoupled  Riccati  equations.  Unlike  the 
H 2  methodology  however,  the  optimization  problem  yields  a 
set  of  Riccati  equations  whose  solution  is  stabilizing  only  if 
IT  H  11.  The  reasons  for  this  characteristic  are  very 
complex,  and  outside  the  scope  of  this  paper.  For  a  detailed 
explanation  of  this  fact,  see  Doyle,  et.  al.  13).  This 
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can  be  circumvented 


,  however,  oy  introducing  a 


scaling  factor  y. 

The  y_  Iteration.  The  level  to  which  the  oo-norm  can  be 
reduced  is  limited.  In  other  words,  it  cannot  be  reduced  to 
any  arbitrary  size  desired.  Therefore,  there  exists  some 
optimal  co-norm  obtainable  for  any  given  system,  no  matter  what 
controller  is  used.  Since  the  Riccati  equations  are  designed 
to  return  a  compensator  which  reduces  ||T2W1I00  to  1/  it  makes 
sense  that  this  would  not  be  possible  if  the  optimal  value  of 
II T z  II oo  is  <3reate^  than  1.  In  this  case,  the  set  of  Riccati 
equations  has  no  stabilizing  solution.  Additionally,  it  is 
equally  possible  that  the  optimal  value  of  ijTzwi00  is  less  than 
1,  in  which  case  the  compensator  obtained  via  the  Riccati 
equations  is  not  the  optimal.  Thus,  a  scaling  factor  r  is 
introduced.  By  a  basic  property  of  any  norm,  in  particular 
the  oo-norm. 


I  r  Tzu  |B»  >-  ITZX  (7 

where  r  is  an  arbitrary  constant.  Thus,  to  obtain  a 
compensator  which  achieves  the  optimal  jTzw|0O/  it  is  necessary 
to  divide  the  input  or  the  output  or  a  combination  of  the  two 
by  a  factor  equal  to  the  optimal  ||Tzw80O*  T^e  optimal  value  of 

IT  |  .  however,  is  not  known  beforehand.  For  this  reason,  it 

>  zw'm 

is  necessary  to  use  an  iterative  approach. 

To  search  for  the  optimal  compensator,  an  artificial 
scaling  factor  y  must  be  added,  producing  a  scaled  input  w  and 
scaled  output  z,  such  that 
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w 


=  <rr  -  l  t  .  -j , 

z  =  47"  z  (7.5) 

The  block  diagram  of  the  resultant  system  is  shown  in  Figure 
7.2.  Note  that  the  problem  has  not  really  been  changed,  only 
the  way  in  which  it  is  perceived. 


Figre  7.2.  System  with  Scaled  Inputs  and  Outputs 


Now  solving  the  Riccati  equations  for  the  system  with  w  as  the 
input  and  z  as  the  output  will  produce  a  compensator  such  that 
II T z w  1  oo  is  e<Iual  to  1  •  The  iteration  comes  about  from  the  fact 
that  if  r  is  picked  too  large,  a  less  than  optimal  solution 
will  be  obtained,  and  if  r  is  too  small,  no  stabilizing 
solution  will  be  obtained.  Thus,  it  is  possible  to  determine 
upper  and  lower  bounds  on  the  optimal  |TzwUco  which  may  be 
narrowed  as  much  as  desired,  but  the  optimal  may  never  be 
reached  exactly. 

In  addition  to  the  necessity  to  perform  this  iteration  on 
r,  there  are  other  problem  requirements  and  scalings  needed 
before  the  theory  can  be  applied. 

Enable*  Requirements  and.  Scaling.  To  solve  the  H 

00 

problem  the  following  conditions  must  be  satisfied; 

1)  Dn  =  0  (7.6) 

21  D22  =  0  (7.71 
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(7  .  5  i 


3)  a)  D12l  D12  --  I 

b)  °2 1  D21T  =  1  (7.9) 

4)  (A,B^)  controllable  &  (C^,A)  observable 

5)  (A,B2)  stabilizable  &  (C2,A)  detectable 

Note  that  these  conditions  are  very  similar  to  the  conditions 
for  the  H2  case,  with  the  exception  of  Condition  4,  which  is  a 
stronger  requirement  than  in  the  H2  problem.  Also,  both 
Conditions  1  and  2  may  be  removed  when  necessary  [5]  &  [11], 
Removing  Condition  1,  however,  complicates  the  problem,  and 
Condition  2  is  trivially  satisfied  for  the  problem  studied 
here.  As  described  in  Chapter  VI,  the  modified  system  used  is 
controllable  and  observable,  so  no  additional  modifications 
are  necessary  to  meet  this  requirement.  Condition  1  is 
satisfied  in  Chapter  VI  by  introducing  weights  on  the  inputs 
and  outputs.  Exactly  the  same  weights  as  used  in  Chapter  VI 
are  used  here,  so  no  further  comment  on  these  weightings  is 
warranted . 

As  in  the  H2  problem.  Conditions  3a  and  3b  require  that 
the  input  u  and  the  output  y  be  scaled.  This  scaling  is 
accomplished  exactly  as  before,  via  a  Cholesky  decomposition, 
which  yields  unique  upper  triangular  matrices  that  satisfy 


(7.10) 


(7.11) 


The  scaling  is  appended  to  the  original  plant  G  as  before. 


7-5 


which  together  with  the  y  scaling,  results  in  a  p  1  .a n t  whicn 

looks  as  shown  in  Figure  7.3. 


Figure  7.3.  System  with  Scaled  Inputs  and  Outputs 


With  these  scalings  all  in  place,  the  new  system  becomes 


x=Ax+Blw+B2u 

<v  A/ 

z  =  Cl  x  +  D12  u 
y  =  C2  x  +  D21  w 

where 


(7.12) 

(7.13) 

(7.14) 


(7.15) 

(7.16) 

(7.17) 

(7.18) 

(7.19) 

(7.20) 


Note  that  Dll  and  D22  have  been  excluded,  as  they  are 
identically  zero  to  meet  Conditions  1  and  2  above.  The 
problem  is  now  in  the  form  required  for  application  of  the 
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Solution.  This  method  requires  obtaining  the  real, 
unique,  symmetric  solutions  X  and  Y  to  two  uncoupled  Riccati 

oo  ao 

equations 


where 


F,  X  +  X  F.  -  X  G,  X  +  H  =  0 
loo  ool  oo  1  oo  1 


F*  Y  +  Y  F-  -  Y  G~  Y  +  H0  =  0 
2  ao  oo  2  ao  2  oo  2 


F,  = 


4  -  B2  D12  C1 


°1  =  B2  B2T  -  B1  B1T 


H.  = 


-  T  ~ 
C  C 


C,  = 


(  1  “  °12  °12  }  C1 


F_  = 


A  -  B,  D21*  C2 


**  m  **  ~  m  ~ 

G2  "  C2  C2  -  C1  C1 


H-  = 


B1  B1 


Bi  (  1  -  d21  d21) 


(7.21) 

(7.22) 

(7.23) 

(7.24) 

(7.25) 

(7.26) 

(7.27) 

(7.28) 

(7.29) 

(7.30) 


Note  that  the  above  Riccati  equations  are  not  the  same  as  in 
the  H2  problem.  In  addition  to  the  extra  scaling  on  and 
C1,  the  quadratic  terms  G  (with  appropriate  subscripts)  have 
an  extra  term  in  them.  This  extra  term  is  the  limiting  factor 
which  makes  the  problem  Inherently  more  difficult  than  the 
H2  problem.  For  more  detail,  see  (3)  and  [4]. 

The  optimal  controller  for  the  problem  is  then  given  by 


B, 


+  Y 


oo 


(  c. 


(7.33  ) 


K , 


(7.32) 
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cp 

D 


cp 

where 


C 


0 


v  7  •  ^  i  } 

(7.34) 


Kc  =  (  B2 


K,  =  Y  C,' 
f  oo  2 


There  are  three 
the  solutions  to 
controllers  for 


~  m  *+ 

X  +  D.  -  C. 

oo  12  1 


)  (  I 


21 


T 


additional  conditions  which  must  be 
the  above  Riccati  equations  to  be 
the  system: 


(7.35) 

(7.36) 

met  for 
stabilizing 


1) 

X  > 

CO  ~ 

0 

(7.37) 

2) 

* 

8 

|v 

0 

(7.38) 

3) 

P(X 

00 

Y  )  i  1 

00 

(7.39) 

where  p  denotes  the  spectral  radius,  defined  as 


*<*„V  =■  “*  (ab3  (^/2  <*„*„>])  (7'401 

where  X  denotes  eigenvalue.  Equation  7.39  guarantees  that  the 

inverse  of  (I  -  X  Y  )  exists,  as  required  for  equation  (7.35) 

oo  oo 

above . 


Elflblea  Set-up 

As  mentioned  previously,  the  identical  problem  set-up 
used  for  the  H2  problem  will  be  used  for  the  problem. 

PcQcedmc 

The  solution  was  obtained  by  taking  the  problem  set-up 

as  described  in  Chapter  VI  and  iterating  values  of  y  until  an 

acceptable  solution  was  obtained,  r  was  reduced  to  a  value  of 

500  which  yielded  solutions  X  and  Y_  which  satisfied  the 

00  00 
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conditions  .above .  The  spectral  radius  of  the  solution 
obtained  Is  0.9254.  Although  somewhat  less  than  the  optimal 
value  of  1,  this  solution  is  very  close  to  optimal,  since 
testing  y  =  499  yields  a  spectral  radius  greater  than  1  (the 
spectral  radius  increases  exponentially  as  the  optimal  y  is 
approached)  and  the  solution  obtained  for  y  =  501  differed 
only  slightly  from  the  solution  for  y  =  500  (the  1-norm  of  the 
difference  in  the  compensators  obtained  was  less  than  10). 

Although  the  same  problem  set-up  was  used  for  the  and 

H  solutions  Dresented  here,  one  problem  was  encountered 
00 

during  the  design  process  which  must  be  mentioned.  The  author 
originally  sought  to  use  bandpass  filters  on  the  disturbances 
inputs  to  more  accurately  reflect  the  expected  noises.  This 
worked  well  for  the  H 2  problem,  but  caused  difficulty  for  the 
problem.  When  attempting  to  obtain  solutions  to  equations 
(7.21)  and  (7.22),  an  error  message  was  generated  which 
indicated  that  the  system  presented  was  uncontrollable  when  y 
was  reduced  below  values  which  yielded  spectral  radii  of 
approximately  0.3.  This  (uncontrollability)  was  known  not  to 
be  the  case  from  prior  analysis,  so  it  was  suspected  that  the 
problem  was  arising  from  some  ill-conditioned  terms  in  the 
problem  posed.  In  order  to  implement  bandpass  filters  and 
still  satisfy  the  conditions  reflected  in  equations  (7.8)  and 
(7.9),  it  was  necessary  to  place  zeros  in  the  filter  which 
were  widely  spaced  in  frequency.  This  is  demonstrated  in 
Figure  7.5,  which  displays  the  Bode  plot  for  one  of  the 
original  filters.  Examination  of  this  figure  reveals  that  in 
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obtain  a  two  decade  passband,  it  is  necessary  to  place  two 
zeros  six  decades  apart,  for  a  difference  in  magnitude  of  106. 
This  large  separation  of  zeros  was  suspected  to  be  the  cause 
of  the  numerical  problems,  and  when  the  bandpass  filters  were 
converted  to  high-pass  filters,  the  problem  completely 
disappeared.  For  consistency  between  solutions,  the  highpass 
filters  were  also  used  for  the  problem,  even  though  the 
bandpass  filters  caused  no  problem  there. 

Order  Reduction  ai.  the  Compensator 

As  in  the  H_  case,  the  H  solution  yields  a  compensator 

Z  oo 

which  is  of  the  same  order  as  the  modified  plant  presented  to 
the  Riccati  equations,  so  model  order  reduction  was  required 
for  the  compensator  obtained  for  this  case  as  well.  This  was 
accomplished  in  a  manner  identical  to  that  presented  in 
Chapter  VI.  A  full  description  of  both  the  full  and  the 
reduced  order  compensators  is  contained  in  Appendix  C. 
Closed-Loop  syatea  RsBsesenta&lQii- 

The  closed-loop  representation  for  the  problem  is 
obtained  the  same  manner  as  described  in  Chapter  VI. 
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VIII.  Per  toruidnce  Curnpar  isons 

Missile  Flygut  Simulation 

To  evaluate  the  performance  of  the  different  controllers 
developed  herein,  it  was  necessary  to  develop  a  digital 
simulation  of  the  missile,  simulating  the  control  system,  the 
changing  environment,  the  noise  inputs,  and  the  target  and 
missile  geometry  calculations.  A  PC-MATLAB  script  file  was 
created  to  carry  out  the  simulation. 

The  simulation  computes  target  position  as  a  function  of 
time,  calls  a  subroutine  which  generates  disturbances  as 
functions  of  time,  computes  the  parameter  variations,  computes 
the  closed-loop  representation  of  the  system  with  compensator, 
computes  the  new  guidance  inputs,  and  increments  the  states 
and  bearing  to  target  for  each  time  increment.  After  end-game 
conditions  are  indicated,  the  simulation  calculates  the  point 
of  closest  approach  and  stores  the  state,  output,  missile 
position,  and  target  position  time  histories  and  the  miss 
distance.  The  simulation  uses  a  0.02  second  time  step,  and 
uses  the  PC-MATLAB  function  c2d  to  convert  the  state-space 
representation  to  a  discrete  iteration  model. 

Ely aut.  Results 

Each  controller  was  evaluated  over  112  flyouts,  composed 
of  16  scenarios  each  of  which  was  run  7  times  under  different 
conditions.  The  7  conditions  were:  one  flyout  in  the  nominal 
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disturbance  inputs  (?,3),  three  flyouts  with  both  parameter 
variations  and  disturbance  inputs  (4,5,6),  and  one  flyout  in 
which  the  target  position  data  is  lost  for  two  0.5  second 
intervals  (7).  The  16  engagement  scenarios  are  described  in 
Table  8.1. 


Table  8.1.  Engagement  Scenarios 


Scenar io 

Velocity 

Trajectory 

1 

200  m/s 

90°  offset 

2 

400  m/s 

90*  offset 

3 

600  m/s 

90*  offset 

4 

800  m/s 

90*  offset 

5 

200  m/s 

45*  offset 

6 

400  m/s 

45*  offset 

7 

600  m/s 

45*  offset 

8 

800  m/s 

45*  offset 

9 

200  m/s 

6g  pull-in 

10 

400  m/s 

6g  pull-in 

11 

600  m/s 

6g  pull-in 

12 

800  m/s 

6g  pull-in 

13 

200  m/s 

6g  pull-out 

14 

400  m/s 

6g  pull-out 

15 

600  m/s 

6g  pull-out 

16 

800  m/ s 

6g  pull-out 
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3 .1-3 .4,  respectively.  The  resultant  miss  distances  (to  the 

nearest  meter)  are  tabulated  in  Tables  3. 2-8. 4,  for  the 

Classical,  H-,  and  H  controllers,  respectively. 

Z  oo 


Table  8.2.  Classical  Controller  Hiss  Distances  (m) 

Environment/System  Condition 

2 

3 

4 

5 

6 

7 

1 

56 

36 

17 

11 

1 

33 

48 

2 

32 

61 

61 

65 

70 

7 

84 

3 

6 

62 

51 

66 

70 

60 

81 

4 

mm 

33 

El 

34 

34 

80 

46 

5 

70 

22 

10 

29 

7 

35 

98 

6 

140 

45 

21 

32 

400 

7 

355 

169 

127 

131 

BE 

BE 

170 

8 

82 

32 

49 

17 

85 

114 

160 

9 

66 

81 

65 

53 

60 

66 

BE 

10 

■a 

45 

El 

55 

50 

El 

IS 

11 

— 

3 

mm 

14 

6 

41 

55 

42 

12 

■a 

16 

2 

8 

IB 

49 

51 

13 

90 

26 

10 

14 

14 

38 

153 

14 

118 

31 

mm 

36 

13 

82 

BE 

15 

102 

7 

31 

2 

El 

54 

8 

16 

57 

20 

17 

24 

24 

3 

3 

Mean 

78 

44 

_ 

35 

37 

79 

96 

Examination  of  the  miss  distance  results  reveals  several 
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Scenario' 12  Engagement 


Table  3.3.  H  ^  Controller  Miss  Distances  (m) 


Environment/System  Condition 

Scenario 

1 

2 

3 

4 

5 

6 

7 

1 

54 

8 

21 

5 

2 

26 

59 

2 

31 

K9 

68 

34 

35 

24 

24 

3 

2 

mm 

42 

42 

42 

22 

m 

4 

16 

9 

11 

mm 

9 

38 

5 

S 

62 

23 

30 

n 

3 

50  i 

123 

6 

EB 

95 

78 

116 

8 

mm 

137 

7 

355 

234 

325 

128 

111 

16 

145 

8 

365 

347 

112 

114 

143 

91 

156 

9 

52 

67 

42 

67 

45 

56 

51 

10 

13 

8 

4 

15 

m 

mm 

17 

11 

24 

12 

12 

24 

ii 

48 

75 

12 

48 

34 

52 

38 

63 

21 

92 

13 

85 

7 

14 

8 

38 

83 

151 

14 

104 

20 

19 

18 

10 

| 

168 

15 

80 

28 

10 

13 

16 

61 

16 

16 

22 

18 

16 

16 

17 

2 

29 

Mean 

91 

62 

54 

42 

35 

64 

82 

surprising  trends.  The  first  is  that  adding  the  parameter 
variations  to  the  simulation  (Conditions  2  and  3)  had  a  very 
beneficial  effect  on  miss  distance  as  compared  to  the  nominal 
case  (Condition  1).  Analysis  of  the  time  history  and 


equations  for  the  parameter  variations  reveals  that  this  vas 
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caused  by  the  fact  that  the  missile  begins  its  flight  in  a 


stable  configuration  when  the  parameter  variations  are 
present.  The  early  trajectory  for  this  condition  exhibits 
less  oscillatory  behavior,  and  thus  achieves  the  required 
heading  correction  earlier  and  with  less  oscillation  than  in 
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tiie  nominal  case.  This  favorable  effect  of  parameter 
variations  was  not  anticipated,  and  complicates  evaluation  of 
the  robustness  of  the  designs. 

Another  surprising  trend  evident  in  the  tables  is  the 
tendency  for  miss  distance  to  not  be  dramatically  affected  by 
noise  inputs.  Figures  8.5  and  8.6  display  the  time  responses 
for  the  H2  controller  in  Scenario  2,  for  both  the  nominal  run 
and  the  run  with  noise  added,  respectively.  The  effect  of  the 
power  plant  noise  is  clearly  evident  in  the  plot,  but  miss 
distance  remains  virtually  unchanged.  This  would  seem  to 
indicate  that  the  controllers  are  more  robust  than  is 
practically  possible,  i.e.  the  more  disturbances  and  parameter 
variations  that  are  added,  the  better  the  system  performs.  It 
is  the  author's  opinion  that  the  parameter  variations 
resulting  from  the  flight  conditions  specified  and  tested 
above  happen  to  fall,  by  coincidence,  in  a  region  advantageous 
to  missile  performance.  This  was  unfortunate  in  this  case 
since  the  purpose  of  the  study  was  to  examine  the  detrimental 
effects  of  these  parameter  variations.  A  more  statistically 
significant  sample  of  parameter  variations  would  have  to  be 
run  in  order  to  truly  get  some  measure  of  system  robustness  to 
parameter  variations  and  disturbances.  Also,  the  time  step 
size  (0.02  seconds)  was  initially  set  based  on  the  nominal 
runs,  where  it  was  maximized  in  order  to  decrease  run  time. 

In  retrospect,  it  is  the  author's  opinion  that  this  time 
step  size  does  not  allow  the  full  effects  of  the  noise  inputs 
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Figure  8.5.  H9  Controller  System  Response,  Nominal  Environment 
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Figure  8.6.  H~  Controller  8ystem  Response,  Noisy  Environment 


to  be  felt,  since  the  PC-MATLAB  function  c2d  implements  a 
zero-order  hold  over  the  time  step. 

Much  the  same  can  be  said  for  the  loss  of  target  runs,  in 
which  it  was  hoped  the  would  prove  to  be  robust  to 
deterministic  disturbances.  Examination  of  the  miss  distance 
data  does  not  reveal  any  significant  advantages  for  the 
controller.  Examination  of  the  state  time  histories,  however, 
does  reveal  some  encouraging  trends.  The  controller  does 
indeed  seem  to  be  handling  the  disturbance  better,  as  the 
transients  are  less  pronounced.  Figure  8.7  shows  the  state 
time  histories  for  scenario  1  using  the  classical  controller, 
and  Figure  8.8  shows  the  state  time  histories  for  the  same 
scenario  using  the  controller .  Notice  that  both  the  body 
rate  and  the  antenna  rate  display  less  severe  transients  for 

the  system  with  the  H  controller. 

00 

While  these  results  do  not  quite  exemplify  the  advantages  of 
designing  controllers  using  the  H 2  and  methodologies, 
sufficient  data  is  available  to  reach  some  qualitative  conclusions 
and  point  the  way  to  future  investigations.  These  will  be 
discussed  in  Chapter  IX. 
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Figure  9.7.  Classical  Controller  Syste.  Pesponse,  Loss  of  Target 


fin  deflection 


Figure  8.8.  H  Controller  Syste*  Response,  Loss  of  Target 


■a  nd  Cj. 


The  design  study  failed  to  reach  the  goal  of  providing 
quantitative  data  which  exemplify  the  robustness  of 
controllers  designed  using  the  and  design  techniques. 
However,  the  data  obtained  and  the  design  process  itself  allow 
several  observations  to  be  made.  In  addition,  several  lessons 
were  learned  which  will  be  of  use  in  any  follow-on  effort  and 
which  will  be  of  interest  to  anyone  working  with  these 
methodologies . 


Rftb.ust.ness  aJL  U.2  and  CQntrQllers 

The  H2  optimal  control  theory  does  indeed  produce 
controllers  which  can  handle  parameter  variations  well,  and 
the  theory  produces  controllers  which  perform  well  also. 
Though  the  quantitative  results  obtained  were  inconclusive, 
the  qualitative  results  described  in  Chapter  VIII  do  display 
some  encouraging  trends.  Neither  the  H2  nor  the  controller 
were  found  to  have  any  shortcomings  in  the  scenarios 
simulated.  In  fact,  both  performed  quite  well.  This  is 
indicative  of  a  sound  approach,  and  a  more  refined  simulation 
and  a  broader  spectrum  of  scenarios  would  probably  yield  more 
meaningful  results.  It  is  the  author  s  opinion  chat  the 
parameter  variation  results  are  realistic,  but  the  simulation 
requires  a  shorter  iteration  time  step  to  fully  realize  the 
effects  of  noise  on  the  system.  This  refinement  i3  suggested 
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an  area  a  r  farther  study. 

Problem  Complexity 

The  complexity  in  setting  up  the  problem  and  executing 

the  design  effort  for  the  and  H  designs  creates  a 

i  ao 

considerable  work  load  for  the  designer,  even  on  a  problem  as 
simple  as  this  one.  It  is  hoped  that  the  methodology 
developed  herein  for  executing  this  process  will  prove  to  be 
one  of  the  most  beneficial  results  of  this  study.  In 
addition,  the  reader  may  benefit  from  some  of  the  lessons 
learned  during  this  design  effort. 

Lessons  Learned 

During  execution  of  the  design  process,  a  problem  was 

encountered  concerning  the  weighting  of  the  disturbance 
inputs.  This  problem  proved  to  be  very  difficult  to  isolate, 
and  was  finally  traced  to  the  use  of  bandpass  filters  to  shape 
these  inputs.  The  problem  was  finally  isolated  by  repeating 
the  entire  design  process,  using  the  design  algorithm. 

By  doing  this,  the  author  was  able  to  isolate  the  problem  as 
occurring  when  the  bandpass  filters  were  introduced.  A 
bandpass  filter  with  a  one  decade  passband  and  two  decade 
rolloffs  at  both  the  low-  and  high-frequency  ends  requires 
zeros  which  differ  by  five  orders  o £  magnitude.  Apparently, 
this  magnitude  difference  causes  ill-conditioning  of  the 
3tate-space  matrices,  which  results  in  a  system  which  looks 
uncontrollable  and  unobservable  to  the  computer  algorithms 
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to  high-pass  solved  the  problem.  The  designer  is  therefore 

cautioned  to  guard  against  this  occurrence. 

Another  lesson  was  learned  concerning  ’ overraod i f ying 1  the 
filters.  In  one  instance,  the  author  was  able  to  get 
perfectly  acceptable  performance  from  the  system  by 
suppressing  the  gain  on  the  body  rate  input  command  to  a  levei 
of  -60  dB  at  low  frequency  and  placing  a  pole  at  0.000001 
rad/sec.  This  configuration  was  extremely  fragile  with 
respect  to  varying  the  other  filter  weights,  however,  as 
system  performance  degraded  rapidly  when  other  weights  were 
modified  to  'tune'  the  performance.  The  bottom  line  on  this 
lesson  Is  not  to  pursue  modifications  that  look  strange,  even 
if  they  seem  at  first  to  yield  good  results. 

Finally,  a  lesson  was  learned  concerning  converting  the 
system  to  the  form  required  for  implementation  of  the  existing 
theory.  During  the  initial  design  effort,  the  original  system 
was  modified  to  meet  the  necessary  conditions  by  using  the 
boresight  error  as  an  input,  rather  than  a  state.  This 
modification  was  required  because  only  body  rate,  normal 
acceleration,  antenna  rate,  and  antenna  position  were  being 
interpreted  as  commanded  values,  so  missile  heading  had  to  be 
excluded  from  the  state  vector  to  preserve  observability. 

This  resulted  in  omission  of  the  boresight  error  resolver 
dynamics,  and  exclusion  of  the  natural  loop  closure  provided 
by  the  resolver.  The  exclusion  of  this  loop  resulted  quite 
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naturally  In  the  stability  problems  caused  by  any  lVop 
closure,  so  that  it  was  possible  to  obtain  compensators  which 
proved  to  be  unstable  when  the  complete  system  dynamics  were 
simulated.  This  problem  was  overcome  by  introducing  a 
fictitious  boresight  error  command  to  the  problem.  By  adding 
this  fictitious  command,  the  full  system  dynamics  could  be 
included,  and  the  stability  problem  disappeared.  Thus,  the 
designer  should  always  attempt  to  add  components  to  the  system 
representation  rather  than  omit  any  of  the  system  dynamics 
just  to  satisfy  the  conditions. 

Areas  F.mth&t  Study 

There  are  many  options  available  for  carrying  out  further 
study  of  this  design  effort.  The  most  promising  of  these  is 
refinement  of  the  simulation.  It  is  suggested  that  the  system 
simulation  be  converted  to  a  Fortran  routine,  so  that  the 
limitations  of  PC-MATLAB  be  circumvented.  Such  a  simulation 
could  implement  a  much  smaller  time  step,  and  thus  bring  out 
the  full  effects  of  noisy  disturbance  inputs.  Also,  such  a 
simulation  could  be  extended  to  five  degrees  of  freedom,  so 
that  motion  in  both  the  pitch  and  yaw  planes  as  well  as  the 
effects  of  gravity  and  drag  could  be  fully  effected. 

Another  area  of  further  study  related  to  the  above  is  the 
simulation  of  a  broader  spectrum  of  flyout  scenarios.  The 
limited  class  of  flyouts  depicted  herein  could  easily  be 
extended  to  a  more  representative  sample  of  possible  target 
encounters.  This  task  would  also  require  conversion  to  a 
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carter  Fortran  furmat  for  the  simulation,  as  each  flyout  in 
the  PC-MATLAB  subroutine  takes  in  excess  of  20  minutes,  due  to 
memory-crunching  and  array  size  requirements. 

Another  area  of  further  study  is  high-frequency 
unmodelled  dynamics.  The  controller  design  methodology  is 
thought  to  be  robust  to  unmodelled  dynamics,  and  an 
investigation  of  this  would  yield  insight  into  the  accuracy  of 
this  conjecture. 

Lastly,  there  are  some  techniques  which  allow  the 
designer  to  tailor  the  design  to  specific  parameter  variations 
expected.  Use  of  these  algorithms  should  improve  robustness 
to  parameter  variations,  as  the  resultant  design  should  be 
less  conservative. 
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Appendix  A-  Weighting  Filters 


The  weighting  filters  derived  for  the  optiaal  controller 
formulations  are  listed  below.  The  filters  are  given  in 
state-space  form  and  as  Bode  magnitude  plots. 


State.  Space  Representation 

Using  the  nomenclature  of  Chapter  VI,  the  weighting 
functions  W^,  W2,  W^,  and  Wz  can  be  expressed  in  state  space 
form  as  follows,  where  the  A,  B,  C,  and  D  subscripts  match  the 
respective  subcripts  on  W. 
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Bode  Magnitude  Plots 

To  provide  the  reader  with  a  graphical  interpretation  of 

the  above  filter  representations/  Figures  A.l  and  A. 2  show  the 

Bode  magnitude  plots  for  the  filtering  of  inputs  1  and  5 

(9  and  9  )  in  W.  (signals  2  thru  4  are  equally  weighted 

mcom  scom 

across  the  entire  spectrum,  so  no  graphical  representation  is 
deemed  necessary).  Similarly,  Figures  A.3-A.6  provide  the 
Bode  magnitude  plots  for  filter  Figures  A.7-A.10  Bode 

magnitude  plots  for  filter  and  A.11-A.15  Bode  magnitude 

plots  for  Wz.  Each  of  these  plots  is  titled  in  reference  to 
the  respective  input,  disturbance,  or  error  filtered. 
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Figure  A. 3.  Accelerometer  Noise  Input  Filter 
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Antenna  Position  Potentiometer  Input  Filter 
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Appendix  B.-  Simulation  Flyout  Time  Histories 

In  this  appendix,  a  representative  sample  of  the  flyout 
time  histories  is  given  in  the  form  of  missile  and  target 
trajectories  and  state  history  diagrams.  The  large  number  of 
flyout  runs  precludes  inclusion  of  all  the  flyout  histories, 
but  this  representative  sample  should  give  the  reader 
sufficient  insight  into  actual  flyout  conditions. 

tthaP-LLe  and  Target  Trajectories 

Figure  B.l  shows  a  typical  missile  and  target  trajectory 
plot  for  scenario  4  (see  Chapter  VIII)  in  which  the  missile  is 
initially  approaching  the  target  broadside  and  the  target  is 
traveling  at  800  m/s. 

Figure  B.2  shows  the  trajectory  plots  for  scenario  8,  in 
which  the  missile  is  initially  approaching  the  target  at  a  45 
degree  offset  from  head-on,  and  the  target  is  traveling  at  800 
m/s  . 

Figure  B.3  shows  the  trajectory  plots  for  scenario  12, 
where  the  missile  is  initially  approaching  the  target  head-on, 
and  the  target  executes  a  6-g  pull-out  away  from  the  missile. 
The  target  is  travelling  at  800  m/s. 

Figure  B.4  shows  the  trajectory  plot  for  scenario  16, 
where  the  missile  is  initially  approaching  the  target 
broadside,  and  the  target  executes  a  6-g  pull-in  toward  the 
missile.  The  target  is  travelling  at  800  m/s. 
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Figure  B.l.  Scenario  4  Missile  and  Target  Trajectories 
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Figure  B.2.  Scenario  8  Missile  and  Target  Trajectories 
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Iii  this  section,  a  representative  sample  of  the  missile 

system  state  time  histories  will  be  given.  All  the  state  time 

histories  given  herein  are  for  simulation  runs  using  the 

classical  controller.  The  state  time  histories  for  the  and 

H  controller  simulation  runs  display  similar  characteristics, 
no 

and  are  nearly  indistinguishable  to  the  unaided  eye,  with  the 
obvious  exception  of  the  additional  states  of  the  respective 
controllers.  Therefore,  only  state  time  histories  for  the 
classical  controller  will  be  shown. 

Figures  B.5  and  B.6  show  the  state  time  histories  for 
scenario  4,  as  described  above.  Similarly,  Figures  B.7-B.12 
show  the  state  time  histories  for  scenarios  8,  12,  and  16,  as 
labe lied  . 
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Figure  B.5.  Airframe  State  Time  History,  Scenario 
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Figure  B.7.  Airframe  State' Time  History,  Scenario 
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Figure  B.8.  Target  Tracker  Time  History,  Scenario 
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Target  Tracker  Time  History,  Scenario  16 


The  state  matrices  for  the  full-order  H2  compensator  are 
given  below.  The  A  matrices  for  the  compensators  listed 
herein  are  very  large,  due  to  the  large  number  of  states 
present,  and  so  must  be  represented  as  partitioned  matrices. 
Also,  scientific  notation  had  to  be  used  due  to  the  many  large 
and  small  magnitudes  of  the  numbers  contained.  This  makes  the 
matrices  difficult  to  comprehend  at  a  glance,  but  no  other 
representation  is  possible.  The  full  order  compensator  can 
be  represented  by  the  state-space  description 


where  y  and  u  correspond  to  the  plant  error  vector  and  control 
P  P 

input  vector,  respectively,  as  described  in  Chapter  6,  and  the 


other  terms  are  given  by 
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where  the  above  are  given  by 
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5 . 70e-08 
3 . 29e-09 
-6 . 47e-07 
-7 . 10e-06 
2 . 92e+00 
1 . 58e-06 
-4 . 90e-06 
1 . 64e-07 
8 . 53e-02 
1 . 66e-01 
0 

1 . 08e-14 
-1.06e-12 
0 

1 . 17e-16 
-3 . 00e+02 
0 


(C.ll) 


1 . 07e-06 
2 . 54e-07 
1 . 35e-04 
5 . 37e-08 
6 . 79e-06 
1 . 89e-06 
1 . 36e-04 
3 . 20e-09 
9 . 99e-03 
1 . 15e-08 
8 . 51e-07 
1 . 94e-06 
8 . 40e+00 
4 . 12e-07 
3 . 89e-08 
3 . 12e-03 
1 . 64e-03 
0 

6 . 08e-14 
5. Ole-13 
0 

9 . 16e-16 
2 . 04e-13 
4 . Q0e+01 


(C.12) 


C2 


cp2f 


C3 


cp2  £ 


C4 


cp2  f 


C5 


cp2£ 


■[ 


-6 , 9  4 e - 0 1 

2 . 1 3  e  +  j  1 

2  .  3  3  e  +  0  0 

6  .  5  7  e  -  0  i 

-  5 . 0  3  e  -  3  4 

1 . 3  oe -  0  3 

- o . b  7e -  0  7 

-4  .  14e  -0  4 

-9  .  35e  -  0  9 

3  .  4  o  e  -  0  7 

(  C  .  1 3  ) 

-3 . 74e-10 

-1 . 55e -0  5 

5 . 7  5e  +  01 

8 . 39e  +  02 

1.61e+01 

-2 . 26e+01 

-5 . 92e+03 

-  3 . 50e-02 

4 . 73e  +  06 

4 . 08e-05 

( C .  1 4  ) 

-5 . 6  3e-02 

8  .  48e-09 

-5 . 40e-04 

-2  .  lle-03 

6  .  37e-01 

-2 . 67e-09 

8 . 14e  +  0 1 

-4 . 06e  +  03 

1. 43e-06 

-5 . 48e-09 

(  C  .  1 5  ) 

-2.52e-13 

-2 . 9  5e-0  3 

-  3 . 7  5e  +  0  2 

1 . 35e  +  00 

1 . 1 4e -0 1 

-1 . 05e-02 

-1 . 03e-02 

-  3 . 9 le -0  4 

-5.92e-Q8 

4  . 7  6e-10 

( C . 16  ) 

2 . 67e-01 

1 . 2  3e -07 

5 . 01e-08 

1 . 9 le-04 

( C . 17 ) 

-1. 82e-04 

-6 . 12e+02 

-6 . 81e+01 

4 . 07e  +  03 

D1 


cp2£ 


■[ 


0 

0 


0 

0 


0  0 
0  0 


(0.18) 


a.2  R&dussd.  Order  Compensator 

The  above  compensator  can  be  reduced  in  size  by  use  of  the 
Hankel  singular  value  elimination  technique  described  in  Chapter 
VI.  The  Hankel  singular  values  for  the  H2  compensator  are 
given  in  Table  C.l.  The  last  13  Hankel  singular  values  were 
eliminated  to  produce  the  reduced  order  compensator,  as 
indicated  by  the  dashed  line  in  Table  C.l.  The 
resultant  reduced  order  compensator,  composed  of  11  states, 
can  be  represented  by 


xcp2r 

u 

P 


Acp2r  xcp2r  +  Bcp2r  ^p 
Ccp2r  Xcp2r  +  Dcp2r  yp 


( C .  19  ) 

( C. 20 ) 


C-5 


T a )j  1  c  C ,  1 .  H 2  COuipCiiS-a t Of  Hankcl  3  1  nyUlai  v.^Iijr-^, 


4 . 0  2e  +  01 
3 . 15e+01 
9 . 00e+00 
3 . 0  8e  +  0  0 
2 . 69e+00 
2 . 0  6e  +  0  0 
1.98e+00 
5 . 56e-01 
4 . 56e-Ql 
3.54e-01 
1.25e-01 

8 . 76e -0  2 
4 . 27e-02 
2 . 53e  =  0  2 
1. 25e-02 
7 . 30e-03 
2 . 26e-03 
9 . 07e-04 
3 . 43e-04 
1 . 69e-04 
7 . 03e-06 
3 . 71e-06 
3 . 41e-07 
1 . 88e-08 


where 


xcp2r 

Acp2r 

B  o 
cp2r 

Ccp2r 

Dcp2r 


compensator  states 


Islcp2r'A2 
1  B1.p2r  1 
1C1cp2rlc2 


1  D1cp2r  1 


cp2r 


cp2r 


] 


] 


( C .  21 ) 
(  C .  2  2  ) 
(C.  23  ) 
(C.  24  ) 


C-6 


and 


A1 


cp2r 


A2 


cp2r 


A3 


cp2r 


B1 


cp2r 


2 . 61e+05 
1.52e+05 
-2 . 20e  +  0  3 
5 . 18e+Q4 
-3 . 16e+03 
l.lle+04 
-2. 41e+04 
8 . 15e  +  0  3 
-3 . 14e+03 
8 . 6  3e  +  0  2 
-3 . 70e+03 


1 . 14e+05 
6 . 69e+Q4 
-9 . 70e+02 
2 . 27e+04 
-1. 38e+03 
4 . 87e+03 
-1 . 0  5e  +  04 
3 . 57e+03 
-1 . 37e+03 
3 . 75e+02 
-1 . 62e+03 


-7 . 19e+04 
-4 . 19e+04 
6 . 07e+02 
-1. 42e+04 
8 . 69e+02 
-3 . 06e+03 
6 . 6  4e  +  0  3 
-2 . 24e+03 
8 . 66e+02 
-2 . 37e+02 
1. 02e+03 


1. 15e+0S 
6 . 73e+04 
-9 . 84e+02 
2 . 29e+04 
-1.38e+03 
4 . 85e+03 
-1.06e+04 
3. 59e+03 
-1. 38e403 
3 . 73e+02 
-1 . 63e+03 


5 . 71e+04 
3 . 33e+04 
-5 . 18e+02 
1 . 13e  +  0  4 
-6 . 31e+02 
2.173+03 
-5 . 15e+03 
1.786403 
-7 . 07e+02 
1.436402 
-8 . 13e+02 

-4 . 75e-03 
7 . 37e-02 
-1.526400 
-9.95e-02 
4 . 73e+00 
1.556400 
1 . 54e+00 
1.266-02 
4.02er00 
-5. 46e+00 
-2.77e-02 


-1.16e+05 
-6 . 79e+04 
9 . 66e+Q2 
-2 . 31e+04 
1.44e+03 
-5. 08e+03 
1 . 07e+04 
-3 . 63e+03 
1.396403 
-4 . 02e+02 
1.65e403 


8.94e+03 
5 . 21e+03 
1.82e+01 
1.716403 
-2 . 66e+02 
9 . 53e+02 
-1.086403 
2 . 78e+02 
-6 . 49e+01 
1.016402 
-1.216+02 

-1.39e-01 
-1.02e-01 
9.146-01 
-1.196400 
-3 . 56e+00 
9 . 12e+00 
-5 . 07e+00 
-2 . 7  2e-02 
-5.21e-01 
1.34e+00 
1.08e-01 


-1 . 4  4e  +  05 
-8 . 41e+04 
1 . 23e+03 
-2 . 86e+Q4 
1. 72e+03 
-6 . 03e+03 
1. 32e+04 
-4. 50e+03 
1 . 7  3e  +  03 
-4 . 57e+02 
2 . 04e+03 


6 . 07e+06 
3 . 54e+06 
-5 . 13e+04 
1.20e406 
-7 . 35e+04 
2 . 58e+05 
-5 . 59e+05 
1.896405 
-7 . 28e+04 
2 . 00e+04 
-8 . 61e+04 

-4 . 40e-02 
-2 . 78e-02 
3 . 40e-02 
-1.20e-01 
4 . 78e-Q3 
8 . 6  5e-01 
-3.00e-01 
9 . 44e-04 
3 . 15e-01 
7 . 08e-01 
5. 18e-03 


-1.03e+04 
-6 . 01e+03 
8 . 9  4e  +  01 
-2 . 04e+03 
1.20e+02 
-4 . 22e+02 
9 . 42e+02 
-3 . 21e+02 
1 . 24e+02 
-3 . 13e+01 
1. 46e+02 


3 . 19e+03 
1.866403 
-2 . 68e+01 
6 . 34e+02 
-3 . 83e+01 
1.356402 
-2 . 9  5e  +  02 
-9 . 96e+01 
-3 . 77e+01 
1 . 05e+01 
-4 . 57e+01 


9 . 43e+00 
5 . 24e+00 
-1.16e-02 
1.746400 
-1.59e-02 
2 . 97e-01 
-1.166400 
-9 . 23e+01 
2 . 50e-01 
6 . 79e-02 
-3.13e-01 


(  C .  25  ) 


(  C .  26  ) 


(  C .  27  ) 


2 . 01e+01 
3 . 53e+01 
2 . 48e+00 
5 . 03e-01 
6 . 89e-01 
1.136400 
6 . 27e+00 
2 . 92e-01 
1.726400 
1.45e-01 
2 . 97e-01 

(C . 28 ) 


C-7 


Cl 


C2 


C3 


Dl 


3 . 81e-02 

-4  .  68e-02 

-7  .  45e-0i 

9 . 04e-0 1 

(C. 29 ) 

cp2r 

1 . 07e  +  03 

-  3  .  9  le  +  0  3 

-1 . 71e  +  03 

1 .  54e  +  02 

-  3 . 7 le  +  OQ 

-6 . 01e  +  00 

5 . 58e  +  00  - 

-2.78e-02 

( C . 30  ) 

cp2r 

-1 . 72e+03 

1.74e+03 

2 . 16e  +  03  - 

•4 . 7  8e  +  01 

'cp2r 

’  -1 . 47e+01 
-8 . 56e  +  02 

3 . 88e  +  01 
-1. 34e+02 

-1.58e+01 
-9 . 09e  +  04 

(C. 31 ) 

cp2r 

'  0  0  0 

0  0  0 

0  0  1 

0  0  J 

( C . 32 ) 

Full,  Order  Compensator 

The  full  order  compensator  developed  for  the  H  design  can  be 

00 

represented  in  state-space  form  as 

Xcpif  =  Acplf  Xcpif  +  Bcpif  Yp  (C. 33 ) 

up  "  Ccpif  Xcpif  +  Dcplf  Yp  (C.34) 

where  yp  and  up  correspond  to  the  plant  error  vector  and 
control  input  vector  as  before,  and  the  other  terms  are  given 
by 


coif 

cpl  f 

=  compensator  states 

=  [A1cpif|A2cplfIA3cpifIA4cplfIA5cpif] 

( C .  35  ) 

cpi  f 

=  [  B1cpif  1 

(C. 36  ) 

cpif 

(C1cpif  |C2cplf  lC3cpif  ,C4cpif  ,,,;:5cpif  1 

(C.  37 ) 

cpif 

-  1  D1cplf  1 

( C . 38  ) 

where 


C-8 


cpi  f 


'cpif 


-1 .  30e  +  02 
-3  .  32e  +  00 
-3 . 55e-01 
-9 . 61e-01 
3 . 51e-01 
-3 . 92e-02 
1.00e-06 
-2 . 85e-02 
2 . 73e-09 
-3 . 21e-01 
1.13e+01 
-3. 21e-06 
2 . 68e-07 
■1. 33e-01 
3 . 50e -01 
-7 . 49e-06 
1 . 88e-01 
0 

8 . 24e-13 
-1 . 61e-ll 
0 

8 . 60e-16 
■1 . 54e-13 
0 


-2 . 60e-02 
-6 . 86e-08 
1 . 61e-06 
2  .  lle-07 
•8 . 78e-04 
■2 • 75e+0 3 
1.00e+00 
5 . 70e-08 
3. 29e-09 
■6 . 47e-07 
■7 . 10e-06 
2 . 9  2e  +  00 
1 . 58e-06 
•4 . 90e-06 
1 . 64e-07 
8 . 53e-02 
1 . 66e-01 
0 

•4  -  85e-14 
■5 . 25e-12 
0 

2 . 04e-16 
•5 . 68e-14 
0 


1 . 69e+03 
-5. 46e-01 
8 . 02e-01 
5. 49e-01 
1 . 97e-01 
1 . 87e  +  00 
1 . 88e-Q6 
7 . 78e-02 
2 . 03e-09 
6 . 07e  +  00 
-3 . 50e-01 
7 . 94e-Q7 
4 . 18e-08 
-6 . 15e-02 
4 .  G2e-01 
-4 . 00e-06 
-5 . 77e-03 
0 

-3 . 12e-ll 
7 . 39e-10 
0 

-2. 34e-14 
4 . 60e-12 
0 


-2 . 91e+02 
2 . 54e-07 
1 . 35e-04 
5 . 37e-08 
-2 . 07e  +  01 
■S . 67e+Q6 
■1 . 36e-04 
•3 . 20e-09 
9 . 99e-03 
1 . 15e-08 
■8 . 51e-07 
-1 . 9  4e-06 
8 . 40e  +  00 
4 . 12e-07 
■3 . 89e-08 
3 . 12e-03 
1 . 64e-03 
0 

8 . 08e-13 
•4 . 10e-ll 
0 

•4 . 65e-15 
■1 . 73e-12 
0 


4  .  01e+02 
0 
0 
0 

■2. 07e+01 
1.17e+03 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


2 . 22e+04 
5 . 46e  +  00 
1.97e+00 
4 . 50e+00 
-1 . 97e+00 
■1.10e+05 
-1 . 88e-05 
■8 . 78e-01 
■2 . 03e-08 
■6 . 07e+01 
3 . 50e  +  00 
■7 . 9  4e-06 
■4 . 18e-07 
6 . 15e-01 
■4 . 02e  +  00 
4 . 00e-05 
5 . 77e-02 
0 

•3 . 12e-10 
■7. 39e-09 
0 

■2 . 34e-13 
■4 . 60e-li 
0 


5 . 06e  +  0 1 
-2 . 69e-01 
-7 .  lOe-Oi 
- 1 . 9  2e  +  00 
7 . 02e-01 
2 . 58e-02 
2 . 01e-06 
-5 . 70e-02 
5 . 47e-09 
-6. 42e-01 
2 . 27e  +  01 
-6. 43e-06 
5. 37e-07 
-2 . S5e-01 
7 . 01e-01 
-1 . 49e-05 
3 . 76e-01 
0 

1 . 64e-12 
■3 . 22e-ll 
0 

1 . 72e-15 
■3 . 09e-13 
0 


3.21c+05 
-2 . 03e-04 
-1 . 08e-01 
-4 . 30e-05 
5.43C-03 
5 . 34e  +  09 
1 . 09e-01 
2 . 56e-06 
-8 . 00e+00 
-9 . 23e-06 
6 . 80e-04 
1.55C-03 
-6.72e+03 
-3 . 29e-04 
3 .  lle-05 
-2. 49c+00 
■1 . 31e+00 
0 

■S . 46e-10 
3 . 28c-08 
0 

-3.72e-12 
1 . 38e-09 
0 


-4.84e-02 
4 . 07e-01 
7 . 26e-01 
2 . 50e-01 
-1 . 20e+01 
-3 . 44e-01 
-6 . 05e-06 
3 . 29e-02 
1. 35e-09 
1 . 20e-01 
2 . 03e-01 
-1 . 98e-05 
3. 41e-07 
-2 . 14e  +  00 
1 . 05e-02 
-1 . 06e-03 
4 . 71e-01 
0 

-2 . 05e-12 
3 . 85e-ll 
0 

1. 36e-14 
-2.77e-12 
0 

( C . 39  ) 

1 . 22e  +  03 
-5 . 46e-01 
-1.97e-01 
-4.50e-01 
1 . 97e-01 
-1 . 16e+02 
1 . 88e-06 
7 . 78e-02 
?.03e-09 
-3 . 92e+00 
-3 . 50e-01 
7 . 94e-07 
4 . 18e-08 
-6 . 15e-02 
4 . 02e-01 
-4 . 00e-06 
-5 . 77e-03 
0 

-3 . 12e-ll 
-7 . 39e-10 
0 

-2 . 34e-14 
-4 . 60e-12 
0 

( C . 40  ) 


C-9 


A3 


cp  i  £ 


A4 


cpi  f 


-4 . 3  3e  +  00 

1 .  2 5e  +  00 

-2 . 57e+02 

-  1 . 9  3e  -  0 1 

4  .  9 Oe  +  0 1 

-5. 89e-01 

-6 . 8  6e -  0  8 

2  .  54e-Q7 

0 

1 . OOe+OQ 

-3 . 55e-01 

1 . 61e-06 

1 .  35e-04 

0 

0 

-4 . 51e-01 

2. lle-07 

5 . 37e-08 

0 

1.00e+00 

3 . 51e-01 

-8 . 78e-04 

-6.79e-06 

2 . 07e  +  01 

0 

-6 . 63e-03 

2 . 7  3e  +  0  4 

-5.8  8e  +  0  6 

-1 . 42e  +  00 

-1. 53e-02 

1.00e-06 

-8 . 62e-07 

-1. 36e-04 

0 

0 

-2 . 85e-02 

5. 70e-08 

-3 . 20e-09 

0 

0 

2 . 73e-09 

3 . 29e-09 

9 . 99e-03 

0 

0 

-3 . 21e-01 

-6. 47e-07 

1 . 15e-08 

0 

0 

-2 . 86e  +  01 

-7 . 10e-06 

-8. 51e-07 

0 

0 

-3 . 21e-06 

-7.07e+00 

-1 . 9  4e-06 

0 

0 

2 . 68e-07 

1. 58e-C6 

-i . 15e+01 

0 

0 

-1. 33e-01 

-4 . 90e-06 

4 . 12e-07 

-5 . OOe+OQ 

0 

3 . 50e-01 

1 . 64e-07 

-3 . 89e-08 

0 

-1 .  C0e  +  00 

-7. 49e-06 

8 . 53e-02 

3 . 12e-03 

0 

0 

1. 88e-01 

1 . 66e -01 

1 . 64e-03 

0 

0 

0 

0 

0 

0 

0 

8 . 24e-13 

-4 . 85e-14 

8 . 08e-13 

0 

0 

-1.61e-ll 

-5 . 25e-12 

-4 . 10e-ll 

0 

0 

0 

0 

0 

0 

0 

8 . 60e-16 

2 . 04e-16 

-4 . 65e-15 

0 

0 

-1 . 54e-13 

-5 . 68e-14 

-1 . 7  3e-12 

0 

0 

0 

0 

0 

0 

0 

(C.41) 

-1 . 20e-03 

7 . 73e-01 

-2 . 87e+04 

1 . 04e  +  02 

8 . 85e  +  00 

0 

0 

0 

-9 . 28e-21 

-9 . 28e-18 

0 

0 

0 

4 . 83e-22 

4 . 83e-19 

0 

0 

0 

-2 . 22e-21 

-2 . 22e-18 

0 

1.00e+00 

0 

7 . 64e-21 

7 . 64e  -18 

-5 . 13e-02 

-5. 22e-02 

-l.lle+03 

-1 . 65e-01 

-1 . 33e-03 

0 

0 

0 

2 . 56e-21 

2 . 56e-18 

0 

0 

0 

-2 . 72e-22 

-2 . 72e-19 

0 

0 

0 

-3 . 48e-23 

-3 . 48e-20 

0 

0 

0 

1 . 77e-21 

1 . 77e-18 

0 

0 

0 

7 . 00e-20 

7 . 00e-17 

0 

0 

0 

7 . 79e-21 

7 . 79e-18 

0 

0 

0 

7 . 63e-21 

7 . 63e-18 

0 

0 

0 

-8.30e-21 

-8 . 30e-18 

0 

0 

0 

8 . 50e-22 

8 . 50e-19 

-5 . 00e+01 

0 

0 

2 . 20e-20 

2.20e-17 

0 

-2 .00e  +  02 

0 

-7 .15C-20 

-7 . 15e-17 

0 

0 

- 1 . 00e  +  00 

0 

0 

0 

0 

0 

2 . 55e-22 

1.00e+00 

0 

0 

0 

-3 . 00e+03 

-i.30e+02 

0 

0 

0 

0 

0 

0 

0 

0 

7 . 77e- 25 

7 . 77e-22 

0 

0 

0 

-1 . 55e-22 

-1 . 55e-19 

0 

0 

0 

0 

0 

(C.  42  ) 


C-10 


1 

2 . 44e+01 

-  3 . 4  3e  +00 

7 . 8  5e -0 1 

2 . 7  Oe  +  0  2 

0 

-1. Ole-15 

-1. Ole-14 

0 

1 

0 

-4 . 07e-16 

-4 . Q7e-15 

0 

■ 

0 

1 . 45e-16 

1. 45e-15 

0 

0 

1 . 12e-14 

1 . 12e-13 

0 

■ 

1 . 80ef 02 

-6 . 9  4e  +  0  4 

1.01e+04 

6 . 25e  +  06 

1 

0 

2 . 76e-17 

2 . 7  6e-16 

0 

0 

-4 . 68e-16 

-4 . 68e-15 

0 

0 

-7 . 44e-18 

-7 . 44e-17 

0 

1 

0 

2 . 21e-15 

2.21e-14 

0 

■ 

0 

1 . 41e-14 

1 . 41e-13 

0 

0 

4 . 17e-16 

4 . 17e-15 

0 

1 

A5cpif 

0 

2 . 55e-15 

2 . 55e-14 

0 

( C  .  4  3  ) 

■ 

0 

-7 . 32e-16 

-7 . 32e-15 

0 

0 

9 . 12e-16 

9 . 12e-15 

0 

■ 

0 

5 . 21e-15 

5 . 21e-14 

0 

1 

0 

3 . 70e-14 

3 . 70e-13 

0 

0 

0 

0 

0 

0 

-9 . 90e-17 

-9 . 90e-16 

0 

1 

0 

-2 . 19e-16 

-2.19e-15 

0 

■ 

1.00e-02 

0 

0 

0 

0 

5. 31e-19 

1.00e+00 

0 

■ 

0 

-2 . 00e+03 

-2 . 10e  +  02 

0 

1 

0 

0 

0 

-4 . 00e  +  01 
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Reduced  Order  Compensator 

As  in  the  case,  the  compensator  obtained  from  the 
design  was  reduced  via  Hankel  singular  value  elimination.  The 
Hankel  singular  values  for  the  full-order  compensator  are 
given  in  Table  C.2.  As  in  the  H2  case,  the  last  13  singular 
values  were  eliminated,  as  indicated  by  the  dashed  line  in 
Table  C.2.  The  resultant  reduced  order  compensator,  composed 
of  11  states,  can  be  represented  by 
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Comparison  Oi.  gall-  and.  Reduced-Order  Compensators 

To  gain  some  insight  into  the  effect  of  using  the  reduced 
order  compensators  above,  it  is  necessary  to  compare  the 
singular  value  plots  of  the  respective  compensators.  Figure 
C.l  shows  the  singular  value  plot  for  the  full-order  H 2 
compensator,  and  Figure  C.2  is  the  singular  value  plot  for  the 
reduced-order  compensator.  Figure  C.3  shows  these  two  plots 
superimposed.  A3  can  be  seen  from  the  plots,  the  only  major 
difference  occurs  at  very  low  frequency,  where  there  are  some 
directions  of  the  input  vector  which  are  magnified  more  for 
the  reduced-order  compensator  than  for  the  full-order  one. 

This  difference  would  generate  some  steady-state  error  if  the 
simulation  were  run  for  a  very  long  time.  The  short  period  of 
the  simulation  (eight  seconds  maximum),  makes  this  effect 
Inconsequential,  so  the  reduced  order  compensator  can  be  used. 

Likewise,  the  singular  value  plots  for  the  full-  and 
reduced-order  compensators,  shown  in  Figures  C.4  and  C.5, 
respectively,  display  a  similar  trend  when  superimposed,  as 
seen  in  Figure  C.6.  As  in  the  Hj  case,  this  difference  is 
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Reduced-Order  H-  Compensator  Singular  Value  Plot 
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Figure  C.3.  Full-  and  Reduced-Order  H_  Compensators  S 
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Figure  C.4.  Full-Order  H  Compensator  Singular  Value  Plot 
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Judged  inconsequent  Lai  for  the  short  period  of  simulation 

herein  used,  and  the  reduced -order  compensator  can  be  used. 
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